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Abstract. Motivated by the theory of superconductivity and more precisely 
by the problem of the onset of superconductivity in dimension two, many pa- 
pers devoted to the analysis in a semi-classical regime of the lowest eigenvalue 
of the Schrodinger operator with magnetic field have appeared recently. Here 
we would like to mention the works by BernofT-Sternbcrg, Lu-Pan, Del Pino- 
Felmer-Sternberg and Helffer-Morame and also Bauman-Phillips-Tang for the 
case of a disc. In the present paper we settle one important part of this ques- 
tion completely by proving an asymptotic expansion to all orders for low-lying 
eigenvalues for generic domains. The word 'generic' means in this context that 
the curvature of the boundary of the domain has a unique non-degenerate max- 
imum. 
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1. Introduction 

The object of study in this paper is a magnetic Schrodinger operator with Neu- 
mann boundary conditions in a smooth, bounded domain fl. We are interested in 
finding an accurate description of the eigenvalues near the bottom of the spectrum. 
In particular, we will improve estimates given in [IIeMo2. in the case of constant 
magnetic field. 

Apart from its intrinsic mathematical interest, this question is important for ap- 
plications to superconductivity. Precise knowledge of the lowest eigenvalues of this 
magnetic Schrodinger operator is crucial for a detailed description of the nucleation 
of superconductivity (on the boundary) for superconductors of Type II and for ac- 
curate estimates of the critical field i?c3- These applications will be the subject 
of further work and will be published elsewhere. We refer the reader to the works 
of Bernoff-Sternberg |He5t] . Lu-Pan [LuPalL ITuFal ILuPaSL ILuPa4| , and Helf fer- 
Pan HePa for further discussion of this subject and to |TiTi| and |S-JSaTh) for 
the physical motivation. 

Let us fix the notations. The domain C is supposed to be smooth, bounded 
and simply connected. Points (xi, X2) in are denoted by z or x. At each point z 
of the boundary, we denote by the interior unit normal vector to the boundary 
of il. We define the magnetic Neumann operator Ti. by 

V{H) 3u^nu = Hh.nu = {-ihS/^ - A{z)fu{z) . (1.1) 

Here A{z) = (—2:2/2, a;i/2), so that curl A = 1, and the domain ViTi) of the 
operator Ti is defined by 

ViU) = {m e H'^{n) I V ■ {~ihV^ - A{z))u\g^ = 0} . 

The case of the half-plane, = M x M+, will be important for fixing notations. 
After a gauge transformation and a partial Fourier transformation we get, in this 
case and with h — 1, the family of models on the half- line: 

= + (.X + , (1.2) 

on i^(R+) and with Neumann boundary conditions at a; = 0. Important results 
about the operators H^'^ will be recalled in Appendix IXI here we only define the 
notation that will be used throughout the text. Let /t*-^''(C) be the lowest eigenvalue 
of H^ i. Then ^ ^ A^^HO has a unique minimum Op attained at a point that we 
will denote by ^q. The corresponding unique positive, normalized eigenfunction of 
j{N,io denoted by uq. We also introduce : 

^ ^ (1.3) 

The main result of the paper gives the asymptotic expansion of the lowest eigen- 
values of Ti. We define fj,^'^''{h) to be the n-th eigenvalue of Ti., in particular, 

,(1 



and prove the following result. 
Theorem 1.1. 

Suppose that fl is a .smooth bounded domain, that its curvature dfl 3 s t-^ k(s) at 
the boundary has a unique maximum, 

k{s) < k(so) =: fc,„ax , for all s ^ Sq , (1.4) 
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and that the maximum is non-degenerate, i.e. 

k2 ■■= -k"(so) ^ . (1.5) 

Then, for all n G N, there exists a sequence {Cj"'}j^i C K (which can be calcu- 
lated recursively to any order) such that //^"•'(/i) admits the following asymptotic 
expansion (for h \ 0) : 

oo 

~ Qoh - + ciej/^yif (2n - + /,i5/8 j2 h='\^;'> . 

3=0 

(1.6) 

Remark 1.2. 

The semiclassical limit h \, is clearly equivalent to a large magnetic field limit, 
since 



-ihV,- BA{z))u{z)\^dz ^ / \{~ij.V , ~ A{z))u{z)\'^ dz 



In 

Remark 1.3. 

Previous results on the bottom of the spectrum ofTih^n were obtained in |HeMo2| . 
who gave the two first terms in the expansion of ^^^\h) (see HcMo2. Theorems 10.3 
and 11.1];.- 

^i^^\h) = e„h - k^^^c^h^'^ + o{h^'^) . (1.7) 

Remark 1.4. 

It is rather reasonable to believe that the proof of Theorem can be adapted for 
getting a similar result under the weaker assumption that there exists J G N, such 
that 

«;(2j)(so) = 0, for j = 1,2,..., J -I, 

K^^'Ks,)^o, ^^-^^ 

i.e. the maximum is non-degenerate of order 2 J. However we will not pursue this 
further. 

If the uniqueness condition in (|1.4|l is replaced by the assumption that there is 
a finite number of maxima (for which 11. 5|) is assumed to hold), we expect the 
existence of sequences of eigenvalues z'^"'\h) corresponding to each maximum. This 
also follows from the techniques applied in the present paper with a little extra work. 

For applications to bifurcations from the normal state in superconductivity it 
seems important to calculate the splitting between the ground and first excited 
states of TL{h). Let us define 

/^{h) = ^l''^\h) - ti^^\h) . (1.9) 

Corollary 1.5. 

Under the hypothesis from Theorem M.lX A(ft,) admits the following asymptotics : 

A{h) ^ Ciey^ Vefcl/i'/^ + h^'/^ £ h^/'^, . (1.10) 

where = cf ^ - C]" • 

The case where $7 is a disc has been analyzed in great detail in |BaPhTaj , using 
the radial symmetry to reduce the problem to ordinary differential equations. In 
this case the splitting A{h) turns out to become zero for a sequence of values of h 
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tending to 0. This is a complication in the analysis of bifurcation. Thus, in some 
sense, the more 'generic' situation considered in this paper has a nicer property. 
We recall that for the disc it is reasonable to conjecture from jBaPhTa) that : 

A(h) , A(h) 
= hmmf , „ < limsup , „ < +oo . 

We recall also that in the case of a domain with a unique corner, with a sufficiently 
small angle, one has f |Bon2j . |BonDaj ) : 

Mh) 
liminf — > . 

h^O h 



In our case, (|1.1U|) implies 



lim^>0 



hi 

Of course (see Bonnaillie |Bonl| for a discussion inspired by Helffer-Sjostrand 
|HeSjl| [HeSj2| ), if there are multiple minima and symmetries, one expects an ex- 
ponentially small gap between the two lowest eigenvalues. 

The plan of the paper is as follows. In Section |21 we prove a simple non-optimal 
upper bound to the ground state energy. This calculation motivates the more 
systematic treatment in Section O where we introduce a 'Grushin problem' in 
order to reduce the analysis to an effective model on the boundary. The effective 
model allows us to construct quasimodes whose energy corresponds to the lowest 
eigenvalues of Ti to any order in h. Thus we get the upper bound inherent in 
Theorem 11.11 In order to prove that the Grushin approach also gives a lower 
bound, we need to prove suitable localization results in phase space. That is carried 
through in sections01and[Sl Finally, in Sectional we finish the proof of Theorem ll.il 
Appendix fXl recalls a number of results from the analysis of the half-plane model 
that are needed in the calculations. Appendix ^ contains definitions concerning the 
coordinate system near the boundary in which all the calculations will take place. 

2. A SIMPLE UPPER BOUND TO THE GROUND STATE ENERGY 

This section contains a simple variational estimate of the ground state energy 
li^^\h). The motivation for giving this result is a number of remarks and calcu- 
lations appearing in the literature. It turns out that the 'obvious' choice of trial 
functions does not give as good energy estimates as one might expect. This moti- 
vates the more systematic approach in later sections. 

Recall that we have defined the constants Qq and Ci in the introduction. 

Theorem 2.1. 

Suppose n is a smooth bounded domain. Let 

fcmax = supk(s) = maxK(s) , 

s s 

be the maximal curvature of the boundary and let 

k2 = inf (-n"{s)) . 

Then the ground state energy of the operator Ti (defined in (|l.l|l ) satisfies 



hmsup/.-^/4{M(^^(/^) - (e,h k„,M'/' + ^l^hy^) } < 
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Remark 2.2. 

Theorem \2. 1\ does not give the correct coefficient to the h'^^'^-term (compare with 
Theorem M.l]) . The trial function used in the proof below is too simple since it only 
uses the ground state uq in the normal variable. This seems to disprove a belief 
stated in a remark in del Pino, Felmer and Sternberg jPiFeSt| . See also Remark \3.S\ 
below. 
Proof. 

The proof consists of an exphcit calculation with a suitably chosen test function. 
(This is the same test function used in the remark in IPiFeStp . 

Let us consider a point xq on the boundary dVL such that the curvature of dO. at 
is fcmax, the maximal curvature of the boundary. We choose our boundary coor- 
dinates (s, t) (see Appendix Ib)| such that xq has coordinates (0, 0). Let x G C^(M) 
be a standard cut-off function : 

X{t) = 1 for \t\ < 1/2 , and suppx C (-1, 1) . 
Consider now the test function 

0^(s, t; h) = Mt, s; h) xi2s/\dn\) x(tAo) , (2.1) 
where, for a > to be chosen below, 

Mt, h) {2a)'^'h-'/''e"'''/''''\'^°'/''''\o{h-'/h) . (2.2) 

and tQ is the constant from Appendix IbI defining the tubular neighborhood of the 
boundary on which one may use boundary coordinates. 

We will get an upper bound to the ground state energy of the Neumann problem 
by calculating the Rayleigh quotient ^"^jl^y/^ for a suitable in the domain of Ti 
Dili.). Actually, one could also work with (j) in the form domain of the corresponding 
quadratic form g-^. 

After a gauge transformation, we can assume that in our boundary coordinates (see 
Appendix : 

A^{A,,A2)^{-t{l-iK{s)),Q^. 

From now on, we fix the gauge such that this property is satisfied. 
Then 

(</>|H0)=/ / \\{hDt-A2)4>\^ + {l-tn{s)r^\{hD,^AM\ 

X (1 — tn{s)) dsdt . 

Now, using the decay properties of uq and the exponential decay of the Gaussian, 
we first get : 

l.\dn\/2 .oo 

{ct^\n^)= / \\{hDt)U^ + {l-tK{s))-\hDs^AM'\ 

J-\di1\/2J0 ^ ' 

X {l-tK{s))x{'2s/\dn\f x{t/U)f dsdt 

+ 0{h°°) . 



6 



S. FOURNAIS AND B. HELFFER 



Again using the properties of uq (see (|A.15ll l and of the Gaussian, we get 



J-\dn\l2 Jo 



-2asyh'^\,f,, -1/2.^,2 



\an\/2 Jo 

\dn\/2 foo 



h\u„{h-''H)\\l~tHi{s))dsdt 



-|ari|/2 Jo 

X {l-tK{s)y^x{'is/\dn\fx{t/tofdsdt 



0{h°°) 



(2.3) 



It is then clear that by interpreting (1 — tK{s)) ^ as J2n>o ^"^(5)" and computing 
term by term, the cut-off function in t does not affect the computation modulo 
0{h°°). So we get 



(0 1 n^) - 



(2.4) 



\dn\/2 poo 

/ e-^''''/'''^\\u'o{h-^/h)f{l~tK{s))dsdt 
-\dn\/2 Jo 

/■\dn\/2 POO 

/i-5/8V2^ / / e-^"'""'"" \uo[h-^/H)\^y. 

J-\dn\/2 Jo 

i^''^io + i2ash^'^ + t{l - iK(s))p(^rK(s)") x{'^s/\dn\f dsdt 



The next step is to replace k{s) by its Taylor expansion k^'^'^{s) at 0, which leads 
to the equahty (modulo 0{h°° j) : 



(0 1 m 



/i-OO pGO 
/ e-^"'" h\u'„{h-'^h)\\l-tK^-yis)) dsdt 
-00 Jo 



-00 Jo 

\2 



X 



/l^/2^0+i2rf/4 + i(l_ i^Tay^^))! ^ t^K^^Y (g)"^ rfg rfi . (2.5) 



Here the cut-off functions have completely disappeared and the integration in the 
s variable is now over (—00, -|-cx)). 

We omit in what follows the reference to Taylor expansions written in superscript 
"Tay" for k and we use for shortness (1 — tK{s))^^ instead of 
next computations. 

With the change of variables a — -\/2a/i~^/*s, r = h^^^^t, we can continue the 
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calculation as 

(0 I - 

+ h ' 



+00 pOO 



00 JQ 



00 Jo 



e-^Vo(r)|^|eo +r(l - t^-^^^'^y {H^^)) + ^V2^h'/^a 



/2a 



X (1 _ dadr 
with (using that k'(0) = 0, since k(0) = /cmax) 

/ + OC> /-OO 
e""' / \u'oiT)f + {^o + Tf\uo{T)fdTda, 
-00 Jo 



+ OC /'OO 



00 JO 

+ 00 pOQ 



2 / 1 J 

- rK(r)p(«:(0) + -«"(0)- 



2a 



T4 = h^/^ 



00 JQ 

+ 00 /'OO 



00 v/o 

+00 poo 



dr da 



e-- |^,o(r)p(Co + r)r2(«(0) + i«"(0)^^) 



dr da 



Therefore, up to 0{h « ), we get the equivalence 

(</) I Hcj)) ~ h{So + + h^/^Sy^} , 

with 



t\u'o{t)\' dT + T{Co+Ty\uoiT)\'dT 



So^Qoj da , 
= n{0) I e-"^ da 



5*3/4 = 2a I a^e da 



4a 



-e-"' da 



\^o+r)\uo{T)fdT 



T\u'oiT)fdT+ I T{^o+Ty\uoiT)\-'dT 



r'(Co+T)|uo(T)pdT 



From the known moments of uq (see Lemma I A . II below or Fournais-Helffer [FoHel 
(6.15), (6.16) and (6.17)]) we have 

1 

2^ 



T|uo(r)pdr = 
t\u'o{t)\^ dr = C\ + 



r(Co+r)2K(r)pdr = ^(Ci+e^/^), 



o3/2 

2 ' 



^^(eo + r)K(r)pdr = ^ + e^/^ 
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So with Iq = J e '^^ da, h = J cr^e '^^ da, we get 
So = ©0-^0 , 

r o"^^'^ 1 

s,/2 = «(o)/o [ - (Ci + -|-) + -{c, + el/') - + e^/^) 
= -K(o)/o(Ci + e^/') , 

Therefore, we finally find 

I n^) = h@oio - h^/^K{o)io{Ci + e^/') 

+ h'/% [2a - ^(Ci + )] + ) . 

We now compute the asymptotics of ||0||2- Along the same lines as the previous 
computations and with the same conventions, we obtain 



\mi 



-OO /'■CO 



-I 



-OO Jo 

+00 

e~ 

-OO 

+ 00 



e-^ |uo(t)|2(1 - h'/\K{!i^)) dadT + Ounif(/i°°) 
(l - h'/^K{^) t\uo{t)\' dr) da + 0„nif(/i°°) 
e-' (l - h'/^K{^)^) da + Ou„if(/i°°) 



/o - /iV2ye^K(0)7o - /z3/4^f!!M72 + o{hi) . 

4a 



So the Rayleigh quotient becomes 



In 



4q: 



12 ^0 
Since the curvature k has a maximum at s = 0, we see that k"{0) < 0. We recall 
that (j) depends on a and that we can now optimize over a. We recover first the 
fact that the term in 0{hX) is obtained without having to specify a. In the case 
when k2 = —k"{Q) ^ 0, which is our main interest, the optimal choice of a is 



a 



k2Ci 



and we get 



I0II2 V 2 /o 



In the case where k"{0) = 0, we can choose a as small as we wish and therefore 



get 



Using 



i^i^ = Boh - K(0)Ci/i3/2 + oih'/^) . 



Io = J e-"^ da = V^ , l2= J a'^e-"^ = ^ 

we therefore get the result of the theorem. 



□ 
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Remark 2.3. 

In the case where k2 = 0, one would expect that the error term a could be 

replaced by (the stronger) 0{k^) for some s G (7/4,2] depending on the order to 
which the Taylor expansion of k(s) — k(0) vanishes at 0. We will not pursue this 
further. See however also Remark \4.6] 

3. Grushin type approach for upper bounds 
3.1. Main statements. 

In this section we will prove the following accurate result. 
Theorem 3.1. 

Let f2 satisfy the assumptions of Theorem M.lV and let n G N. There exist a sequence 

{Cj"''}j^o *- ^ '^^'^ sequence of functions {<t>j^^}JL[) *^ LF'{^) such that, for all 
N > 0, there exists M > such that, if 



3=0 

(3.1) 



M 

bthx,h)^Y.^^/'cl,f\x), (3.2) 



then (for h\ 0) 



II(w-^m)</'mIIl^ = 0(/^^)I!'^mIU- (3.3) 

With the notations of the theorem, we define (h) as the asymptotic sum 



3=0 

(3.4) 



Consequently, z^^j^ (h) is the truncated sum of z!£\h) at rank M. 



Remark 3.2. 

The lowest approximate eigenvalue z^-^^h) agrees with the calculation from Bernoff- 
Sternberg jBeSt| (see also "SP j up to the order that they calculate (term of order 
/i^/^j. However, it is different from the result stated in del Pino, Felmer and Stern- 
berg [PiFeStI Remark 4.2] that they claim to be sharp. 



Since the operator Ti. is self-adjoint, we can deduce the existence of eigenvalues 
ir the points 
Corollary 3.3. 

Lei n e Pi 
such that 



near the points with asymptotics z'£\ 

Let n G N, M e N and let z^j(^\h) be as above. Then there exist C > and ho > 



dist(z^"^(/i),Spec(7^)) < Ch^ , Vh e (0,/io] . 

Proof. 

This is clear by the Spectral Theorem. □ 
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Remark 3.4. 

In particular, the upper hound announced in Theorem is a direct consequence 
of Theorem \3.1\ 

Proof of Theorem \S.l\ 

The proof is fairly long; so we will split it in different steps described in the next 
subsections. From now on we will assume that the maximum of k, fcmax , is attained 
at s = . 

3.2. Expanding operators in fractional powers of h. 

From |HeMo2. (B.8)] we get that the operator Ti. in boundary coordinates becomes 
n = \{hDs - Ai)a-\hD, - ii) + {hDt - A2)a{hDt - ia)] , (3.5) 

with 

a(s,t) = 1 - tK(s) . (3.6) 

Remark 3.5. 

The representation ofTi given in ()3.5|l is only defined on functions with support in 
[0,to) X [— |(9r2|/2, +|9ri|/2]. We will only apply our operator on functions which are 
a product of cut-off functions with functions in the form of linear combination of 
terms like W^wih^is, h~^t), withw miS(]RxIR+). These functions are consequently 
0{h°°) outside a fixed neighborhood of (0,0). This is similar to the calculations in 
the previous section. We will do the computations formally in the sense that : 

• Everything is determined modulo 0{h°°); 

• a^^{s,t) will be replaced by Xln>o(^'*'('^))") 

• k{s) will be replaced by its Taylor's expansion. 

For any n and N, we will find M and construct trial functions (f)^j^^ (expressed 
in boundary coordinates {s,t) and in the form (|3.2(l ^. localized near (s,t) = (0,0) 
and satisfying 

11(7^ - z'm^)&\\l^ = Oih'')\\&h. , {hD, - i2)0i"^|(,,,)=(,,„) = . (3.7) 

By changing back to the original coordinates, this clearly implies (|3.3|l . We will 
omit the tilda's in the following and thus denote by the trial function in boundary 
coordinates. 

As in the preceding section (see also Appendix IB|| . we choose a gauge where 
ii = -to2(s,t), ^2=0; a2{s,t) ^l-tK{s)/2 . 

We make the scaling r = h^^^^t, a — h^^^^s. Then Ti becomes 

P ^ d-\h'^/^D^ + h^/^T&2)a-\h'^/^D^ + h^'^Th2) + ha-^DrCiDr , (3.8) 

with 

a(a, r) = 1 - h^^^TKih^/'^a) , a2(a, r) = 1 - h^^^TKih^/'^a)/2 . (3.9) 

Thus 

h-'^P = h~^{h^'^D„ + Th2)h^^{h^'^D„ + TO2) + a^^DrhDr . 
We now define 
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and get, after removing the tilda's from the a's, 

P = a-' ((r + Co) + /^'/'^. - r(l - a2))a-\{T + Co) + h'^'D, - r(l - as)) 

+ a-^DraDr~Qo . (3.10) 

We assume that k, is C°° and has a non-degenerate maximum at s = 0. Then, in 
the sense of asymptotic series in powers of /is , we obtain 

a{a, r) = 1 - h^/^TK{h^/^a) = 1 - h^'\n{Q) - tY h^/^+i/^^l^-lM , (3.11) 
and 

m 



2 2 y 1 1 ! 1 

J=2 

From the asymptotics of a and 02, we get (remember the definition of fcs from p.SIl ') 

a{a, t)-i = 1 + h'/\>iiO) - Th^/*^ + Oih^^) , 
a(a,T)-2 = 1 + 2h'/^TK{0) - Th^/^a^k2 + 0{h'/'') , 
-r(l - a^ia, r)) = _/ii/2^2^ ^ r^/i^/V^^ + 0(/i^/«) . (3.13) 
Thus, we can write 

P = Po + h^^^Pi + h^''^P2 + /i^/^Pj + h'/^Q(h) , (3.14) 

where 

Po -i?' + (T + eo)'-eo , (3.15) 
Pi = 2i?,(T + eo) , (3.16) 

P2 = -2r2^(r + Co) + 2T«;(0)(r + Co)' + <Q){rDl - D^tD,) 

= K(0)(2T(T + eo)'-T2(r + eo))+i«;(0)i?r , (3.17) 

P^^Dl- ra'Ur + ?o)^ + 2rV2^(r + ^o) - ^(rZ?^ - Z?,ri^,) 

= i?^ - (2r(r + e„)^ - (r + ^o)) ^ - ^-^^Dr , (3.18) 

and where Q{h) admits a complete expansion : 

00 

3=0 

We define 6P by 

^P = P-Po, (3.19) 

We search for functions 0'^") (h) having an asymptotic expansion in h^^^ and such 
that 

(P- + ®"^ )^(")(fe)^0, (/i; a, 0) = 0. (3.20) 
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The constructed functions will have sufficient decay properties to allow interpreting 
(|3.2Q|I in the sense and therefore, after multiplying by the cutoff appearing in 
(EiH), we get (|X7j) (which implies (j^ ). 
3.3. Reduction to the boundary. 

We will now explain the strategy initiated by Grushin [Gruj (and reference therein) 
and Sjostrand in the context of hypoellipticity [^. Here we use this strategy 
for producing good trial functions and thereby results for the magnetic Neumann 
Laplacian. 

Let us define the operators Rq , Rq and Eq by : 

R+ : 5(M,) ^ 5(R, x (l^),) (3.21) 

0(cr) i-> (/)(cr)wo(T) = ® Mo , 
i?o : S{M.a X (17)r) S(Ra) (3.22) 

f ^ fi<^,T)uo{T)dT , 

"'0 

(3.23) 



Here we abused notation and considered Pq as an operator on L ((R+)t-) in order 
to define Eq. That E^ respects the Schwartz space S{Ra x (ffi+)r) follows from 
Lemma IA.5I 

Notice that Rq is the Hilbertian adjoint of Rq (seen as an operator from 
L'^iRo- X R+;dadT) into ^^(M^)). On the other hand (P - z) is for z e R for- 
mally selfadjoint for the original (/i-dependent) scalar product inherited from 
the change of variable z i-^ {s,t) i-^ (c, r) (that is associated to the measure 
(1 — h~^TK{h~^(j)dadT)) but not for the usual associated to the standard 
Lebesgue measure dadr. 

With the above notations we define matrices of operators 

These operators act on S{M.a x (M+)t-) x 5 (Mo-). Actually we should better think of 
operators applied to formal expansions in suitable fractional powers of h with coef- 
ficients in these S spaces. These infinite formal expansions will then be truncated 
at a suitable rank for defining our quasimodes. So we prefer to write formal expan- 
sions to infinite order, having in mind that we could actually go back to truncated 
expansions if we want a given, arbitrarily small, remainder estimate. 
We note first that : 



Rq 



oSq^I. (3.25) 



An easy calculation gives then : 



V{z)£o = I + IC, 
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where 

_ f{SP - z)Eo {SP - z)R+ 
\ 

If we look for z — z{h) satisfying 

z{h) - ^Zf/is , (3.26) 

e>3 

and having in mind the expansion H3.14|l . we observe that {6P — z) = 0{h^^^) , 
when acting on a fixed function in 5(IR.cr x (K"'")r) and can actually be expanded in 
powers of /is, starting from — z^). So, if we define 

0=0 

then the operator is well defined (after reordering) as a formal expansion in powers 
of /is and 

r{z)£oQoo-I. (3.27) 

Now, 

j^j ^ f}{SP - z)EoV [{5P - z)EoV-HSP - z)R+ 
and therefore, if we write 



we get, in the sense of formal expansions in powers of /is , 

{P~ z)Eoo{z) + R^E^{z) ^1 , (3.28a) 

(P - z)E+{z) + R+E^(z) ^ , (3.28b) 

R^Eooiz)^0, (3.28c) 

RoE+{z)^l. (3.28d) 

So, in particular, if ct i— > (/)oo(o'; h) is a function^ such that 

i?±(z)0oo^O, (3.29) 

then inserting (poo in Ij3.28b|l (i.e. inserting ( 5* j in H3.27|l ). we find 

{P-z)E+{z)^^^0, (3.30) 
where everything is well defined modulo 0{h°°) . 



^More precisely ipooi ■ \ h) ^ EjGN ^^'Pji') ' ^° ^^^^ computations have to be expanded in 
powers of h s . 
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3.4. Construction of trial functions. 

From the above, we see that E^{z) is the following asymptotic series, 

oo 

E^i^) = Y.^-iyRo[{sp - z)Eoy-Hsp - z)K ■ (3-31) 

i=i 

We look as before for 



OO 

j=o 

such that 

Et{zooih))^oo{cT;h)^0, (3.32) 

in the sense of asymptotic series in powers of /i « . Here the functions (jjj are supposed 
to be in 5(R„). 

Lemma 3.6. 

For each n £ N, there exists a unique solution to equation 

in the sense of asymptotic series, and such that 



3=0 

Conversely, for any pair (z{h),ip{h)) such that ()3.32|l is satisfied, with 
z(h) ^ Chi + hi J2j>oQ^^ ^"■'^ ^ X]j>o^^'^J' there exist n and 

c{h) ~ Cjhi such that z{h) = z^^^^h) and ip{h) = c{h)(j)(''\h) . 



Proof of Lemma \S.f)\ 
Let us write 

C30 

3=0 

The terms in this sum will be given in (|3.34() , H3.35|l , H3.38() , and (|3.45|) below. Using 
the definitions (|3.14|) and (|3.19|l . and the fact that R^^ Eq — and EqR^ = 0, we 
get modulo terms of order 0{h's ) , 

E^{z^{h)) ~ - zi) + h^/\P^ - Z2) + h^^\P3 - Z3))R^ 

+ hi Bo P1E0P1R+ + 0{h''/^ ) . 
Since also Rq[t + £,o)Rq = , wc find, using I3.16|l . 

^i = -i?o(Pi-zi)i?+ = zi. (3.34) 
Furthermore, using again H3.17|l . 

E2 = -RoiP2 - Z2)Rt =Z2- k(0)(/i,1 + /l,2) , 

with 

h.i = / [2r(r + eo)' - r2(r + ^o)]uI{t) dr , 
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and 

/•OO 

h,2 = i / uo{T)DrUo{T) dr . 

Jo 

Using Proposition IA.2I we get 

^1,1 + Ii,2 — —C\ , 

and therefore 

E2 = Z2 + k(Q)Ci . (3.35) 
The term E3 becomes, inserting Pi and P3 from (|3.1()|) and (|3.18|) . 

= Z3-Dl + ^ih,i + 11,2) + iDlh , (3.36) 
where we have introduced 

l2 = / {T + Co)Mr)Po'{r + ^„)uo{T)dT . (3.37) 
Jo 

Using Proposition IA.2l and Proposition lA.3l we have 

1-4/2 = 3Civ/0^, /i,i +-^1,2 = -Ci , 

and we therefore get 

I, 2 

E3 = Z3'- 3Ci Ve~oDl -Ci^. (3.38) 

Remember that k(s) has a non-degenerate maximum at s = 0, so k2 — — k"(0) > 0. 
The first terms. 

In order to get the equation (|3.32|) to be satisfied, we choose 

zi=0 , Z2 = -k(0)Ci , (3.39) 

which imphes 

Ei^O , E2^0. (3.40) 
With this choice, (|3.32|l becomes 

^h^/^E3<l>o + Oih''^) . 
So we determine 23 and (t)Q by 

Sa'/'o = . (3.41) 
Let us solve the equation E^4) = . It reads, with k2 = — k"(0) , 

— -,2 , fc2cr% , 



So, after the scahng s — y gj^g o" i "^c find that Z3 should be an eigenvalue of the 
harmonic oscillator 



Thus the possible values of Z3 are : 



4"^ = Ci y/Qok2 (2n - 1) , where n e N \ {0} . (3.42) 
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In particular, using the inequality 3Cl^/Qo — 1 ^ 4/2 < 1 (see Proposition lA.Sp . 
we get that Zg^' is smaller than the value in Theorem 12. II 
Remark 3.7. 

A second look at the calculations above (comparing with Section\^ shows why The- 
orem \^n\ does not give the correct ground state energy to order /i^/*. By using a 
trial state which has the simple form (I2.1|l and H2.2|l . we would not see the term 
P\EqPiRq in the first line of and therefore the last term, 4D^/2, would 

he missing in the second line of (|3.3fi(l . Thus the harmonic oscillator discussed 
above would become 

Dl + , 

instead. This harmonic oscillator has ground state energy \J'^^§^ in agreement with 
the result of Theorem \2.1\ 

The iteration procedure. 

Let us define 11 to be the orthogonal projection on {(/)o}^. Notice that 0o depends 
on the n chosen in (|3.42l) even though we do not explicitly recall this dependence 



in the notation. We will choose (t)j such that 

(pj _L (/)o for all j > . (3.43) 

The term of order hi'^i in (j3.32|l becomes 

j-i 

^s-/-, +^i^fe0,-i-fe =0. (3.44) 

Notice that 

F,^Cj~R^QjR+ + F, , (3.45) 

where Fj only depends on ^i, Z2, and {CkYtZ^o- -^Y taking the scalar product with 
in the equation H3.44|l . we therefore get, by using ()3.41fl and the property that 
i?3 is self adjoint, 

O-i||0o|P = (00 .R(,Q,-iRt<t>o) - (00 I ^;--i0o> - (00 I ;^/^fc0,-i-fc) . 

Since 0o 7^ 0, this equation determines Q-i G C as a function of zi^ Z2, z^, {Cfc}i=o 
and {0fc}{Zo • "^^^ property that Q-i G M will be only proved later. 

Upon projecting the equation (|3.44|) on {0o}^5 a-nd remembering the choice 
(jOSji . we get 

UE^ncj^j = -n(^ i^fe0,_i_fc) . (3.46) 

fe=0 

Since nSall is invertible on {0o}^7 H3.46|l together with H3.43|l determines (pj. 
Uniqueness. 

Suppose that 2:1,22 are not given by the choice in H3.39|l . For concreteness, let us 
suppose that zi ^ 0. Then the equation H3.32|l implies that ~ 0. Thus H3.39|l is 
the only nontrivial choice. 
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Furthermore, in the construction above we imposed that _L (j)o for aU j > . 
Suppose we do not impose that condition. Let (/i^ be the solution constructed above 
and let (f>j be the new solution. Then we can write each as 

= <j)'j + Cj(j)a , with (jj'j _L 0o, c-j e C. (3.47) 

We now write 

Hh) ~ '/'o + ^ /i^'^Vj - c{h)^o + J2 '^'^'^i- ' 
with c{h) + ^^'^"^0 ' 

and 0j ± (j)Q . 

By linearity, we therefore find that is the solution (/'o+X!j>i h^^^4'j constructed 
above. 

This finishes the proof of Lemma IT^ □ 

Using Lemmas 13.61 IA.4I and IA.5I we can finish the proof of Theorem 13.11 Let 
be one of the formal solutions from Lemma [3.61 By stopping the 
formal sum at a finite number of terms we obtain partial sums {z'"^\h) , (j}^^\h)) , 
solutions to 

EUz^Z\h))<l>^Mih) = h"RM{h) , (3.48) 

where Efj is also defined by stopping the expansion of E^. Since the 0j's are 
Schwartz functions and all involved operators respect the space S (they are diff'er- 
ential operators whose coefficients are smooth with polynomially bounded deriva- 
tives), the remainder RM{h) in H3.48II is bounded in S. Using Lemma [A. 51 Lem- 
ma and the fact that all terms in E'^j preserve the Schwartz space (differential 

operators with polynomially bounded, smooth derivatives), we sec that E+^(t)^^^{h) 
defines a finite sum, whose coefficients are in the space iS(IR. x R+). Thus, the 
procedure described in Subsection 13.31 above (reduction to the boundary) gives a 
solution ^->M{h) — E'^j(j)^^j {h) to the equation 

{P-ZM{h))i,M{h)^0{h^') . 

Here the right hand side is in S and controlled in 0{h^^) for any semi- norm on 5, 
thus in particular in the norm. 

Moreover, H being selfadjoint, we can now prove that Qj G K. and this finishes 
the proof of Theorem 13.11 □ 



4. Space localization 

In this section we will prove that the ground state is well localized both in s and 
t. In the following Sectional we will prove a similar (slightly weaker) localization 
result in the frequency variable ^ corresponding to s. 
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4.1. Agmon estimates in the normal direction. 

If (/) is a function with compact support in Q, i.e. (f> G C^{Q), then 

/ \{-thW ~ A)cl>\^ dx ^ f \{-thW - A)^\' dx > h\ml,^^,^ = hml,^,,^ . (4.1) 

Since 1 > Go , this implies (compare 1 • h with Oq • h) that functions with energy 
below our upper bounds from Theorem 13.11 cannot be localized in the interior of 
(i.e. away from the boundary), as /i — > . The powerful method of Agmon estimates 
can be applied to strengthen this property into an exponential localization of the 
eigenfunctions (corresponding to the bottom of the spectrum) in a neighborhood 
of the boundary. This is the content of the following proposition. 

Theorem 4.1 (Normal Agmon estimates). 

Let ho > 0, M €z (00,1). Then there exists C,a > and hi g (0, ft-o] such that 
if {uh)h£{o,ho] a family of normalized eigenfunctions ofHh.Q with corresponding 
eigenvalue n{h) satisfying fi{h) < Mh, then, for all h € (0, hi], 

f e^"'^''''^'='^"^/'''^\\uh{x)\^ + h-^\{-ihV - A)uh{x)\^)dx <C . (4.2) 

Proof. 

The proof is similar to (but easier than) the proof of Theorem l4. 91 below. We omit 
the details and refer to Helffer-Morame ^HeMo 2 Section 6.4, p. 621-623] or Hclffcr- 
Pan |HePa| . In |HeMo2| only the ground state is considered, but it is immediate to 
see that the analysis goes through for the eigenfunctions corresponding to higher 
eigenvalues. □ 

As a corollary, we get the weaker but useful estimate for near the boundary. 
Corollary 4.2 (Weak normal Agmon estimates). 

Let the assumptions be as in Theorem \4 1\ For any integer k, there exist C > 
and ho, such that 

\\distix,dn)''uh\\mn) <Ch^ ,yhe (0,/io] . (4.3) 

Remark 4.3. 

The L^ statement in Theorem \4.1\ can be converted to an L°° result using the Sobolev 
imbedding theorem. See HeMo2 Theorem 6.3] for details. 

4.2. First lower bound. 

In order to get good localization properties of the eigenfunctions in the variable 
parallel to the boundary (the s variable), we need to improve the lower bound on 
the ground state energy from (|1.7|l . We will prove the following improvement of 
|HeMo2l Theorem 10.3]. 

Theorem 4.4. Let D, be a bounded region with .smooth boundary satisfying the 
assumptions of Theorem M.lX Then 

f^^'\h) > Qoh - Cik^^^h^'^ + 0{h''/^) . (4.4) 

Proof. Since Vt is bounded, we have /cmax > 0. Using the results from |HeMolj . we 
may localize to the region near boundary points with k{s) > 0, so we may assume 
without loss of generality in the proof that k{s) > for all s. 

The proof of Theorem 14.41 is similar to the proof of HcMol, Theorem 10.3]. 
We just need to do one of the estimates slightly more carefully. In particular, the 
proof goes by comparison with the case of a disc. For a disc, one can calculate the 
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ground state energy with great precision by using the rotational symmetry. This 
was carried through in jBaPhTaj . We state one of their results in the following 
form. 

Theorem 4.5. Let ii^^\h, b, D{0, R)) be the ground state energy of the operator in 
in the case where curl^ = b (independent of x G ft) and fl — D{0,R); the 
disc of radius R. Then there exists C > such that if 

bR^/h > C , 



then 



(1) 



{h, b, D{0, R)) > Qobh ~ Cib^^^h^^^/R - Ch^R-^ . (4.5) 



Notice that in the case of a disc, the curvature k is constant, k — R ^. 

Let p > 0. We can find a sequence {sj,/i}^^Q'' in ]R/|(?r2| and a partition of unity 

{Xj.h}f=o^ on R/\dQ\ such that suppxj^/i n suppxfc.ft = if j ^ {fc - 1, A;, fc + 1} 
(with the convention that N(h) + 1 = 0,0—1 = N{h)). Furthermore, we may 
impose the conditions 

snppxj,k ClSj,H + [-hP,h''], ^x'h = l, El^^^^l'^^^"'"- (4.6) 

j 3 

We will always choose the Sj,h such that \sj^h\ < |i9ri|/2. 
Let xii X2 be a standard partition of unity on R : 

X? + X2 = 1 , supp XI C (-2, 2) , XI = 1 on a nbhd of [-1, 1] . (4.7) 

Let us define 

X3M{s,t) = Xj,h{s)Xi{t/h^) ■ 
We will also consider Xj.h as a function on Q, (by passing to boundary coordinates) 
without changing the notation. For i]j G D{TC), we can write 

i'^T.^lh^ + ^lh^^ (4.8) 

3 

with 

O3AX) = X3{t{x)lh') , for J = 1, 2 . (4.9) 
We get by the 'IMS'-formula : 

I Hi,) = I n{Y.xlh4' + olni^)) (4.10) 

3 

= Y.(X3M^ I nXjMi^) -h" j \Vx,,h?\M^ dx 

3 

+ {e2,h^ I m2,h^) -h^ I iV^s.hPlV'P dx . (4.11) 



In particular, for Tp — u^^"^ being a normalized ground state wave function, we 
get using the weak normal Agmon estimates (in the t variable) and the condition 
p < 1/2 (we will choose p = 1/8 in the end). 

M^'HM = Y.(^3A'^ I -Hx.A'^) + o{h'-^n . (4.12) 
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We will write each of the terms in (• | •) in boundary coordinates and compare with 
the similar term with fixed curvature. Define 

ii(s, t) = -t{l - tK{s)/2) , a(s, t) = 1- tK{s) , 

Bj,h ■■= J e[xj,hu\^^]is,t)dsdt , (4.13) 
e[/] := a-^\{hD, - A,)f\^ + a\hDtf\^ . (4.14) 



and 



with 



Similarly, we define 



and 



with 



Aj.h ■■= J ej,h[Xj,hU}l ]{s,t) ds dt , 



(4.15) 



(4.16) 



We observe that Bj^h — {Xjji^^h'' I '^Xj.h.u^h'') ^nd will compare Bjji and Aj^h- We 
clearly have 



e[Xj,hU^h\s,t) - ejj,[xj.hu^^\s,t) + /i(s,t) + /2(s,i) + h[s,t) , (4.17) 



and 



/3 = 2a-i5j{(ii - A^,,^n)x,.nu'i\hDs - ii,,,;.)x,,/.^ii'^ } • 
Notice that for s G Sj^h + [—hf ^ h^] , we have, since k'{0) = , 

\K{s)-K,,h\ = \s-Sj,h\- f K'{{l-e)s,,h+£s)de <ChP{\s,,h\ + hP) . (4.18) 

Jo 

Thus, 

|a - aj, ,,1 = t\Kis) - K,- ,,1 < ChP{\sj,h\ + hP)t , 

|a~^ - a-}} < ChP{\sj^h\ + hP)t , for t< 2hP . (4.19) 
We estimate, using (|4.19f) . for any e > , 

\fi{s,t)\<ChP{\s,,^\+hP)te[xjWkh{s,t) 

< C'e/ii/4+2''(|s,-, J + e[x,;,.«i'^](s,t) + C'e-i/i-i/4t2e[x,^,.4')](s,t) 

C\h'/^+^P{\sj,h\ + hP)\,AXj^u\^'^]is,t) + C'e-'h-'^H^e[x,,^^^^^^ . 

(4.20) 

(4.21) 



< 



We also estimate /2 and by 

f2{s,t)<Ch^Pt^\x,.^u^^W 
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and 

|/3(s,t)| < 2Ct^hP{\s,^h\ + hP)ajl\xjMui'\hD, - ii,,, 

(4.22) 

Thus, we get by combining H4.17|l with (|4.20l) . 14.21|l . and H4.22|l and integrating, 
B,.h > {1 - Ceh^/^+^P{\s,,h\ + hpy}Aj,h 

- Ce-^h-^/^ j fe[xj,hu'i\s,t) dsdt 

- C(/i2'' + e-i/i-i/4) j t^\x^ hu^^\dsdt . (4.23) 
From Theorem l4.5l we get the estimate 

> (eo/i-Ci(fc„,ax-comin(|s,-,,|2^1))/i3/2_c;j2^||Xj,h4'^ir , (4-24) 

for some cq > , using the non-degeneracy of the maximum. Therefore, using that 
\sj.h\ < |9f2|/2 , we get that, for e sufficiently small and p = 1/8 , 

{1 - Ceh'/''+"'{\s,M\ + hpf}Aj,h > {Ooh - Cifc.„ax/i'/' - Ch'/'')\\xj.hu\Pf. 

(4.25) 

Therefore, we get from (I4.23f) , for the given choice of e and p — 1/8 , 

Y.^^ '- ^ (eo/i- Cifc,„ax/i'/' - (4.26) 



-Ch'^/^ I Yfe[xj,hU^h\s,t) ds dt ^ Ch-^'^ I t^\u^^\^dsdt . 
J . J 

The weak normal Agmon estimates and easy manipulations (as in |HeMo2j ) give 
that the last two terms in (|4.26(l are bounded by Ch'^^^. Therefore, the theorem 
follows from (|4.26|l and H4.12|l , remembering that p = 1/8 . □ 

Remark 4.6. The above proof actually extends to the case where kmax 'is a non- 
degenerate maximum of higher order, i. e. 

Kis) = K,,^ + as^^ + Oi\sr+') , 

with N > 2 and a ^ 0. In that case (|4.18|l becomes 

\k{s) - K,- „| < C/i^ds,- + hP) . (4.27) 

This implies that (|4.23(l becomes 

B,,n >{l- C6;ii/4+2p(|,^. ,|2A^-i + hPf}A,M 

- Ce-^h-^/^ j t^e[xj,hu''yi\s, t) ds dt 

-C(/i2p + £-1/1-1/4) f t^l^^j^u^^^l^dsdt . (4.28) 
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But instead of H4.24|l . we get 

Aj,h > {Qoh-Ci{k^,^-c„max{\s,j,f^,l))h^^^ -Ch^}\\XjMU^,!'^f , (4.29) 

Since |sp^ > js^^^"'^^"'^^ for small s, this implies the result. Actually, in this case 
one should be able to optimize the proof above (in particular choose p < 1/8) and 
get a better error bound than 0{h'^/^) m 

It is convenient to have a lower bound of the operator Ti. in terms of a potential 
Uh- That is our next statement. 

Theorem 4.7. There exist eo , Cq > such that, if 

k{s) /smax - eos^ , (4.30) 

and 

^f t{x)>2h'/\ 
^" \eo/^-Cl«(s)/l3/2-Co/^^/^ if t{x)<2h'/\ ^ ' 

then 

{u\nu)> I Uh{x)\u{x)\'^ dx , (4.32) 
Jo. 

for all u e D{H) and all h G (0, 1] . 

Remark 4.8. It is very likely that one could replace k by k by k in (|4.31() (see also 
[HeMoll Proposition 10.2]J. However, we do not need this improvement. 

Proof With 6i^h , 02,h as in H4.9|l and p = 1/8, we have 

{u I Hu) = {6i,hu I Ti6ijiu) + {62,hU \ 'H92,hu) 

~Ch^/^ f \u\^dx. 

J{/ll/S<t(j;)<2/ll/S} 

Since {62.hu \ H02jiu) > /i||02,/iM|P, it therefore suffices to prove that 

{u\nu)> [ Uhix)\u{x)\^ dx , (4.33) 

for aU u e H^{n) and all h e (0, 1], where 

~ jjh , if t{x) > 2h^/^ , 

^''^''^^\Qoh-C,k{s)h^l^-C'^hV\ if t{x)<2h^/\ ^^-^^^ 

and 7 = (1 + 8o)/2 and Cq is some positive constant. 

Let v}j^^ be a ground state for Ti — Uh with ground state energy fl^^\h). We will 
prove that /l^^) (/i) > . 

Since Go < 7 < 1, the normal Agmon estimates, Theorem 14.11 are also valid for 



Using the 'IMS'-formula, and notations as in the proof of Theorem 14.41 we get 

fi^^Kh) > J2{x,.h4'^ I {n-Uh)x,.hui'^) + I {h-UH)\e2,hu^u''?dx 



Ch^^ \ui'^\'dx 

J{h^/^<t{x)<2h^/^} 
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Modulo choosing C'q sufficiently big, it therefore suffices to prove that 

3 

>-Ch'/^ jY,\x,.Huf^''dx. (4.35) 

j 

Since the normal Agmon estimates hold for u^j^^ , we can now go through the proof 
of Theorem 14.41 with ujl^ replaced everywhere by u]^^ ■ We replace ujl^ everywhere 
by u]^\ in particular in the definition of Aj^h and Bj^h, which are then denoted by 
Aj.h and Bj,h ■ In particular, we get as in 14.2411 

\h > {eoh-Ci{kn,,^-comm{\sj,h\^a))h^^^-Ch^}\\xjMui'''\\^ ■ (4.36) 

So 

>-Chy^XjA'^r ■ (4.37) 

Therefore, 

j 

= E (^^^'^ - i^oh ~ C^k{s)h'^')\\XJJ.u^H^r) 
j 

-Ch-y^ I ^t'e[x,,hU^^\s,t)dsdt-Ch-^/^ I t^\u^^\^dsdt . 
J ^ J 

Using the weak normal Agmon estimates to bound the last terms by 0{h7/^)^ this 
implies H4.35|l and therefore finishes the proof of Theorem 14.71 □ 
4.3. Agmon estimates in the tangential direction. 

Theorem 14.71 can be used to obtain exponential localization estimates in the tan- 
gential (s-)variable. 

Theorem 4.9 (Tangential Agmon estimates). 

Let ho > 0, M > 0. Then there exist C, a > and hi S (O,ft.o]; such that 
if {uh)h£]o,ho] family of normalized eigenf unctions of Ti with corresponding 

eigenvalue ^i{K) satisfying the hound 

fi{h) < Ooh - Cifc„,ax/i^/' + , V/i e (0, ho] , (4.38) 

and if Xi ^ C'o" function from (|4.7() . then, for all h G (0, hi], 

J e''"\'^^^\'^^"\l{t{x)/hy'^){\uh{x)\^ + h-^\{-thV - A{x))uhix)\^} dx < C . 

(4.39) 
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Proof. First we observe that there exists /? > such that, for all S" > 0, we have 
(with X2 from (|4.7|l ) 



Uh 



Ch^ 



V(x2(:^)Xi(t^)) e''^\^^''^'''\uu\'dx 



(4.40) 



But it follows from Theorem 14.71 and (|4.38|) that if a is chosen such that fia^ < 
eoCi/2, then 



(x2(s^)Xi(xw)e 



Uh 



X X2(s)^)xi(7r^)e"l^l'/''''^,. 



> -Co- M)/^^/lx2(-^)xi(T^)e"l^l^/''^^^ 



Uh\ 



(4.41) 



Therefore, it follows from (|4.40|l and H4.41|l that, for a sufficiently small and S 
sufhciently big, 

||X2(^)xi(x^)e"l^l^/''^^^.|r 

<Ch''' I \w{^,{^,)^,{^,))\\^-\^\'/^"'\uu\^ dx . (4.42) 



Now, 



h'" I |v(x2(^)xi(xw))|V"i^i^/''^''Kpdx </ + //, 



(4.43) 



// := C 



/{/!.l/8<t<2/!.l/8} 

On {/ji/s < t < 2h^'^} we have |s| < \dQ.\t/h^'^, and clearly |s| < \dQ.\/2, so we 
have 

II<c( e"l^"l'*/''''V/.|'rf^- 

J{/ii/s<f<2/ii/S} 

By the normal Agmon estimates fTheorem l4.1|l . this implies that 

II <C . 

To estimate /, we use that 1X2(5^178)1^6^"'^' ' is bounded uniformly in h and 
get 



I <C j \uh\'^dx = C 
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Since also Xi( g/Ji/s )e^"l''l^/'''^^^ is bounded uniformly in h, 14.42|l implies that 

||xi(7r^)e"l^l'/''''^.|r<C. (4.44) 

The bound on 

Jn 

in (|4.39() now follows in the same way by inserting 14.44(1 in H4.40|l . □ 
Corollary 4.10 (Weak tangential Agmon estimates). 

Let the assumptions be as in Theorem \4-9\ Let x G C'o(M), suppx C (— ^0,^0); 
with the constant to from the definition of the boundary coordinates in AnvendixT^ 
Then, for all k > 0, there exists C > such that 

f \s{xtx{tix))\uhix)\Ux<Ch''^^ . 
Jn 

The proof of the corollary is immediate. 

5. A PHASE SPACE BOUND 

For our analysis of the low lying eigenvalues of 7i, we need, apart from the 
localizations in s and t, to have a precise localization in Ds- This is the goal of the 
present section. 

5.1. Main statement and main step of the proof. 
Theorem 5.1 (Localization in Dg). 

Let M > 0, /iQ > and xiy X2 S C°°(R) be a standard partition of unity as in (I4.7|l . 
Let (s,<) be the boundary coordinates introduced in Appendix\^ chosen such that 
k(0) = fcmax and let e in (0,3/8). Then for all N > exists Cn > such that if 
i'U'h)hG{o,ho] ^■s a family of normalized eigenfunctions o/7i = 7i(/i) with eigenvalue 
IJ-{h) satisfying 

Kh) < Qoh - Cifc,„ax/i'/' + , (5.1) 

and the operator Ws acting on functions localized near the boundary is defined by 

W.Xi{t/to) = Xi(4g/|g»l)X2( '^'^'^;~^°' )xi(4s/|9»|)xi(4tAo) (5.2) 
with to from l|B.l|l . then 

\\WsXii4^t/to)u4^, < CnH'' . (5.3) 

and 

\{W,xii'^t/to)uh I nih)WsXi{'^t/to)uh) \ < CNh'^ . (5.4) 

Let us be more explicit about the meaning of the operator WsXi{i/U)) ■ On 
the support of i xi(4i/to), we can use boundary coordinates {s,t) (see Ap- 
pendix EJ- Thus, for each 4> £ L'^{Vl), f{s,t) := xi(4t/to)</' is a |9r2|-periodic 
function in s. After multiplication by xi(4s/|9ri|) we find a function with sup- 
port in (— ^■^^) X which we extend by zero to a function (with compact 
support) on K X M+. This function we still denote by xi(4s/|i9r2|)xi(i/io)0- On 

M X IR_|_ the meaning of the operator X2(-' ^/ ' ) is obvious (for example using 

the Fourier transformation). After multiplying a second time by xi(4s/|c^f^|) we 
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get a new function xi{'^s/\dQ\)x2{- f^' )xi(4s/|t^^|)xi(V^o)'/' '^ith support 

in (— , ) X ]R+ which we may reinterpret as a function on a neighborhood of 
the boundary of expressed in boundary coordinates. 

Thus, WsXi(t/to) is an /i-pseudodifferential operator (or rather /i^/^^'^-pseudo- 
differential operator). We will use elementary commutation properties of such op- 
erators. The relevant results (and much more) can be found in introductions to the 
subject, such as |DiSj| and |Ro| . 

Remark 5.2. As a shorter notation, instead of (|5.3(l and (|5.4|l we will write 

\\WsXi{^t/to)uh\\^, + {{WsXii'it/toW I nih)WsXii^t/tn)uh)\ = Ou„if(/i°°) . 

Here the subscript 'unif is included to remind us that the constants (in (|5.3|l and 
(|5.4(l ) are uniform for eigenf unctions in a suitable energy interval ( as given in 

Proof of Theorem \5.1\ 

Let < (5 < 1/2 and define W, Xs, Xs,o, Xt, and xo by 

W := X2( '^'^'^r^°' ) , Xt Xdt/h^^'-') , (5.5) 

Xo:=Xi(4iAo), Xs,o:=Xi(4s/|91^|), x« Xi(s//i'/'^') ■ (5.6) 

We will choose S small such that : 

0<<5< (|-e)/4. (5.7) 

By using the normal Agmon estimates f Theorem 14. Ill , it suffices to prove the fol- 
lowing localized versions of 1)5. 3|l and 1)5. 4|l : 

||VF,Xt«h||L2 -Ounif(/i°°) , (5.8) 

{WsXtUh I n{h)WsXtUh) = O^Mh°") ■ (5.9) 
We start the proof of (|5.8)l and (|5.9(l by the easy identities 

^i{h)\\W,XtUh\\'' = ^{xtW:W,xtUh I nuh) = + T2{uh) , (5.10) 

with^ 

Tiiuh) := {W.xtUh I nWsXtUh) , (5.11) 
T2iuh) := ^{uh I {xtW:WsXtn - 2xtW:nWsXt + HxtW:WsXt)uh) . (5.12) 
We will also use the following estimates 

WxtUhW' < 1 , ixtuh I Hxtuh) < Ch . (5.13) 

Only the second estimate in (|5.13(l deserves comment. It is however an easy conse- 
quence of the standard identity 

Xtnxt = Ux^n + Hx't) + h'\WxtW 

and of the estimate : 

H{h) < Ch , (5.14) 

resulting from Assumption (|5.1|) . 



"^Notice that Tj depends also on a choice of a pair (xt , Ws) and that we will have to consider 
different pairs in the induction argument. 
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Induction argument. 

The proof of H5.8(l and H5.9(l will be obtained by proving by induction that p{N) is 
satisfied for any iV e N, where p{N) is the following statement. 

Statement p{N): 

For any xt and W as in (|5.5(l . then 

WWsXtUkh^ = Ou„if(/i'^(^-^-'^) , (5.15) 
{WsXtUh I 7^(/i)W^.Xt«/.> - Ounii{h^''^^-'-'^+^) ■ (5.16) 
Initialization = 0. 

The estimate (|5.15|) is trivially satisfied for iV = and H5.16|l is a consequence 
of l|5.11|l . H5.1U|I . I|5.14|l . H5.7|l and Proposition 15.71 fProposition 15.71 is somewhat 
stronger than needed at this step). 
From iV to TV + 1. 

Suppose now that we have proved p{N) for some > 0. Given xt and W, choose 
Xt and W satisfying the same assumptions, but being slightly 'larger', i.e. 

XtXt = Xt , ww = w . 

We introduce Ws := Xs,oWxs.o- Then we consider (ph ■= XtWsUh instead of Uh 
and assume p{N), with the pair {Ws,Xt)- We come back to (|5.10|) and observe, 
using the rough ft,-pseudodifferential calculus, that Tj{uh) = Tj{(t)fi) + Ounii{h°°) 
(j — 1,2). We notice also that Proposition 15 . 61 implies that one can take (ph — (p^ 
in Propositions 15.31 and 15.71 Therefore, (|5.1U|) and Proposition 15.71 applied with 
'Ph ~ 4>h together with p{N) leads to (|5.16|l m^i . Finally, Proposition 15.31 and 
(I5.16|l Mj-i give (using (|5.1|) ') ^^^n+i and this finishes the induction. 

Thus, (|5.3|l and H5.4|l are proved and we have reduced the proof of Theorem 15. II 
to the proof of the three Propositions 15.31 15.61 and 15.71 which will be given in the 
next subsections. □ 

5.2. Step 2 : Lovi^er bound for the local energy Ti{Lph)- 

Proposition 5.3. Lei S e C^(M), S = 1 ok [-i, i], S = on M\ [-1, 1]. Suppose 
that e € (0,3/8), that 5 satisfies (15.71) and let C > 0. Then there exists cq > 
(depending also on the constants implicit in 0{h°°) in (|5.19|) and (|5.20|l below) and 
for all N gN exists Cn > such that ififh G D{T-L) satisfies 

supp^ft C {^(a;) < 2/1^/2-5} ^ (5 ^^y) 

{\ifh? + h-\-ihV - A{x))ifh\'') dx<C , (5.18) 

X2is/h'/'-'){\^h\' + h-'\{-zhV - A{x))iph\'} dx = 0{h^) , (5.19) 
and 

4s h^/'^D ~ fn 4s 
Then, 

Tii^h) ■■= {^k I n^k) > (Ooh - Cifc^ax/i'/' + Wvkf - CNh"" . (5.21) 
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Proof of Provosition \5.cl[ The term Ti{Lph) is an integral 

Ti{iph)^ J {a-'\{hD, - Ai)iph\'' + a\ihDt)iph\^) dsdt . (5.22) 

We introduce a localized version of Ti{ifh), 

fiiifh)^ J xl{a-^\ihD, - Ai)iph\^ + a\ihDt)iph\^) dsdt . (5.23) 

Using the localization estimates in s (see (|5.19|) '). we obtain : 

Ti{^h)^fi{^h)+0{h°-) . (5.24) 

On the set 

{|s| < n {<< 2/^1/2-*} , 

we have 

a(s, t) = 1 - tfc„,ax + 0(/i'/*-'^) , a-\s, t) = l+ tkm.. + Oih"^-""^) , 

ii(s,t) --<(!- ifc,nax) + 0(/l'/'"'') . (5.25) 
Therefore, with 

ai(t) l-ifc,nax , a2(i) l + 2tfc„iax , A{t):=-t{l-tk^,^) , (5.26) 
we get, using (|5.24(l . 

Tii^h) > (1 - h^'^-^')Q[^h] + 0(/i'/'*-5')|l(/.„|p + 0{h^) , (5.27) 

where 

j xl{a2{t)\{hD,- A{t))f\^ + \{hDt)f\^)a^{t)dsdt . (5.28) 
It is clear, using again H5.19|l . that we can remove the localization Xs and get 

T,{Vh) > (1 - /.3/^-3*)g[S(^)^,] + 0(/,V4-5.)||5( Ji.)^^||2 ^ (5(/^oo) ^ 

(5.29) 

where 

Q[f]-=£^ [a2{t)\{hDs~A{t))f\^ + \{hDt)f\^)ai{t)dsdt. (5.30) 

Now the coefficients in Q do not depend on s, so we can make a Fourier decompo- 
sition of the quadratic form. Let us define 

f{s,T) ■.^h^/^xi{h'r)f{s,h^/^T) , 
i:=k^^^h^TXi{h^T/2) , 

5c(t) (2^)^1/2 / e-<^/(s, r) . (5.31) 
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Notice that the function i is uniformly bounded on M+. Then 

Qd9]--= {(l + 2/ii/2-^^(r))[/ii/2^ + r(l-/.i/2-^£(r)/2)]'|5(r)|2 

+ y{T)\']{l-h'''-H{T))dT . (5.32) 

The form is the quadratic form on (1 - h^/'^-H) dr) D L'^{R+,t'^{1 - 

^i/2-(5^-j dr) defining a selfadjoint unbounded operator f)(C) on the space (1— 
h^/^~^£)dT) : 

Similarly, we can introduce the quadratic form on iJ^(R+) n L^(]R+, r^dr) : 

q^ig] := {{h'^'C + rf \9{r)\' + \g'{r)f} dr . 

with associated operator t)o(C) on i^(M+,dr) which is the Neumann selfadjoint 
realization of : 

In the two cases, the form domain is the same space and the operator domain 
involves the Neumann condition at r = 0. 

Lemma 5.4. There exists cq , C , M > , such that if\h^/'^C - Col > M/i^/^-^-^/z ^ 
then 

infSpec[)(C) >eo + comin(l,|/ii/2^-eo|'), (5.33) 
and if l/ii/^C - < M/ii/4-3V2^ then 

infSpecf)(C) > {eo + 3Ci|eo|(/i'/'C-Co)'-Cifc„,ax/i'/'} 

-C(|/»i/2c-eoP + /»^/'|/i^/'C-Co|). (5.34) 



Proof of Lemma \5.4\ The proof is similar to a calculation given in |HeMo2l Sec- 
tion 11], so we wiU be rather brief. Since < £ < C, < < /i"* and 5 e (0, 1/2), 
we get for all / in the form domain of f)(C), 

qdf] < (1 + Ch^^^'')q'c[f] + h'/'-''\\fr , (5.35) 

and the same inequality is true by exchanging q^ and q^. Thus, by the variational 
characterization of the eigenvalues, 

\^^,m))-NiUO)\<Ch'^^-''{l + ^^Jii)o{0)} ■ (5.36) 

Here /Xj(f)) denotes the jth eigenvalue of the self-adjoint operator [) (with f) = f)(C) 
or f)o(C)) ■ Now it follows from (|A.4|I that, for some co > , 

m(f)o(C)) > eo + comin(l,|/ii/2C-Co|') , 
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and therefore we get from l|5.36(l , if M is chosen sufhciently big, that 

•^0 ______ i,_1/2a ^ l2\ I7_1/2a I ^ 71 r7_ l/4-3i5/2 



MHO) > + y min (l, \h^/\ - ^oP) for l/i^/^C - Co| > A//i'^ 



(5.37) 



Note that from 1/4 - 35/2 > 0. For \h^^'^( - < Mh^/*'^^/'^ , we will con- 

struct an explicit trial function for [)(C)- With P^^ being the regularized resolvent 
from (|A.14|) we write 

/c(r) = uo{r) - 2{h'/\ ~ i^P^\{t + Co)«oW](t) (5.38) 
+ \[h^l\ - eo)'Po"'{(t + eo)Po"' {{f + ^^)u^{t')\ it) - huo{t)}{T) . 

We note that /.^(t) belongs to the domain of [)(C) and a straightforward calculation 
gives that 



{UO - [eo + {h'^'C - Co)'(i " 4/2)]}/c|| < c\h'^'C ~ ^of 

where I2 is the constant from H3.37|l . i.e. 



(t + $o)uo(r)Po"i [{t + 6)wo(i)](r) dr . 



Define now 



/C := fc - h'/^k„,,.Po' [{—+ 2t(^o + r)^ - r^i^o + t)}uo 
Again a straightforward calculation, using the decomposition 

f)(C) - f)o(C) := hi'' {jrm^.T-r + nh"^C + ^ - h^'^-'lr/if 

gives that 

- [00 + (/i'/'C - eo)'(l - 4/2) + /l^/^fc.„axA^]}/c 

with 



(5.39) 



(5.40) 



(5.41) 



N 



(5.42) 
(5.43) 
(5.44) 



Uo(r){ — + 2r(eo + t)^ - r2(eo + r)}?/o(r) dr . 

Using Propositions IA.2I and IA.31 we get 

N^^Ci, l-4/2 = -3Cifo = 3Ci|fo| . 

This, together with H5.35|l which permits to have a lower bound of fi2{h{()), finishes 
the proof of Lemma [5.41 fsee HeMo2] for a similar argument). We have actually 
obtained the better 

fiiim) - {©0 + 3Ci\^o\{h^/\ - Co)' - Cifc„,ax/j'/'} 

~C{\h'/\-^o\' + h'V^'C-io\). (5.45) 

□ 



Lemma f5 .41 has the following consequence 
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Lemma 5.5. Let (ph,^ md S satisfying the assumptions of Provosition TB^ Then 
there exists cq > such that 

Q[^{^^)^h\>{Q^h-Cik^,^h^'^ + coh^+''') J \E{-^^)iphmi-tk^,^)dsdt 

The proof of Lemma 15.51 is immediate. 
End of the proof of Proposition IsTsl 
Using l|5.17(l and H5.19|l , we get that 



+ 



Therefore, Lemma l5 . 51 imphes that 

g[S(^Vh] > ieoh-C,k^,.h^/^ + coh'+''^Chy*-^')\\^h\\l,^^^ . (5.46) 

Combining (|5.4t)|l with (|5.29|) , and using the choice of 6 from (|5.7|) , yields (|5.21|l . □ 
5.3. Step 3 : Preservation of localization. 

Proposition 5.6. Let e G (0,3/8) and let 5 G (0,^). Then there exist a,C > 
such that if (uh)/ie(o.ho) ^■s ^^^^ family of functions from Theorem \5.1\ and Xti W are 
as in (|5.5|) . then (f)^ := XtWsUh satisfies 

gat(.)AV^||^^|2 ^ ;,-1|(_,;jV _ A{x)M^} dx<C, 
X^(s//ii/'-') I' + h-^\{-ihV - A{x))^r,f} dx = Ounifihn • 



Proof of Provosition \5 .(A 

We only consider the localization in s, since the localization in t is much simpler. 
Let us define 

xl{slh^/^-')\{-ihW - A(x))^u? dx . 

We will only prove that T = Ounif (^°°), the remaining estimate in Proposition 
being easier. We write, with Xs '-^ X2{s/h^^^^^), 

T = Ri+R2, (5.47) 

^i := Jxta~'\{hDs-Ai)^h\^dsdt, R2 := j xla\{hDt)^h? ds dt . 

Since a is a bounded function on suppxt and Wg commutes with Z3(, we have, with 
(• I •) being the inner product on L'^([-\dn\/2, \dn\/2\ x M.-dsdt), 

|i?2| < C{{hDt){xtUh) I W,xlWs{hDt){xtUh)) . (5.48) 
Now, with Xj defined by 

Us)^W^s/h'l^-'), 
WsxlW, = {xl + xl)WsxlW,{xl + xl) ■ (5.49) 
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Since e < 3/8 , we see by repeated commutations that Xi^sxl — Cunif(^°°), so 

m < C{xl(hDt){xtUh) I WsxlW, xl{hDt){xtUh)) 

0^rni{h^)\\{hDt){xtUh)f . (5.50) 

Now the Agmon estimates in s and t easily imply that i?2 — Ounii{h°°). 

The estimate of Ri is similar but slightly more complicated, since Ws does not 
commute with Ai. We estimate 

|i?i| <4'^+4'' , (5-51) 

Ri^ --^C jxl\{hDs+t){xtWsUh)\^dsdt, 

Ri^ -^C J xllAi-tl-'lxtWsUhl^dsdt, 

On suppxt, \Ai — t\ < C, so R\ ' can be controlled like i?2 by Ounifl^"")- 

For r''-^ we use that W commutes with {hDg + 1) and the Agmon estimates in 
s, to write 

= C{{hDs + t)xtuh I WsxlWs{hDs + t)xtUh) + Ou„if(/^°°) ■ (5.52) 

Notice that 

\\{hD, + t)xtUhf < CWihDs + Ai)xtUh\\^ + C < Const . (5.53) 

So we can, like for the control of i?2, localize modulo Ounii{h°°) on the support of 
X2 and get 

< CixlihD, + t)xtUH I WsxlWsxl{hDs + t)xtuh) + Ou„if(/i°°) 

< C'WxlihDs + t)xtuhf + OunAhn 

< C"{\\xl{hDs - A,)xtUh\\^ + WxlxtUkW') + Ou„if(/i°°) 

= Ounitihn ■ (5.54) 

Here we used the tangential Agmon estimates to get the last inequality. This 
finishes the proof of Proposition 15. 61 □ 

5.4. Step 4 : Control of the commutator T2{(ph) • 

5.4.1. Main statement. 

Proposition 5.7. Suppose that e e (0, 3/8) , that S satisfies H5.7|l and let a , C > 

0. Then there exists cq > (depending also on the constants implicit in (|5.5()|) 
below) and for all N Cz N there exists Cn > such that if iph G DiTi.) is such that 

gat(x)//ii/^||^^|2 ^ _ A{x))iph\^} dx < C , (5.55) 

X2is/h'^'-'){\^k\^ + h-'\i-zhV - A{x))iph\'} dx = 0{h^) , (5.56) 
Then, with Wg and Xt from (|5.5|l and 

T2{vh) ■■= li^h I {xtW:WsXtn - 2xtW:HWsXt + nxtW:WsXt)vh) , (5.57) 
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we have 

\T2{vh)\<coh'''^-'-'({^h\n^h)+h\\^h\\'')+CNh'' . (5.58) 

The proof of Proposition is based on successive decompositions of the 'com- 
mutator'. 

5.4.2. First decomposition for T2{iph) ■ 

Since Xt locaUzes near the boundary, we can use boundary coordinates (s, t). Thus, 
we get, with a = 1 — tK{s), 

W: =a-^W,a^W,+W,- t?. x.,oW^Xs,o , W -.^ —[W, k{s)] . (5.59) 

a 

Remember that Ws = XsfiWxsfi- Let S e C^{^) satisfy, S(s)xs,o(s) = Xs,o(s), 
suppS C (— ■^^,-^^). Clearly, T2{iph) = T2{'E,Lph). Now we can calculate the 
'commutator' in T2{^^Ph)', 

xtw:wsxtn - 2xtw:nwsxt + nxtw:wsxt ; 

as an (pseudodifferential) operator on L^(R^), where we extend the curvature func- 
tion k{s) (appearing, for instance, in the expression for Ji in boundary coordinates) 
as a periodic function of s G K. 

Using the localization estimates in s from (|5.56() combined with the fact that xt 
commutes with W and W , we therefore get 

T2{^h) = ^ixloVh I {Ci+C2)xlo^h) + Ounidhn , (5.60) 

Ci:=[xtW,[xtW,n]], C2 :=C2,i+C2,2 + C2,3 , (5.61) 

C2,i:=W[xu[Xun]]W , €2,2:^ Wx^t[W,H], C2.3 := [H, t?] W^Xt • (5.62) 
We will calculate and estimate these commutators. We write 

nih) ^ a-^{hDs - Ai)a-^{hDs ~ M) + a-^{hDt)a{hDt) 

= Hi+H2+H3 + Hi , (5.63) 

with 

Hi := {hDs - Ai)a~\hDs - ii) , 7^2 := -ih^{hD, ~ ii) , 

:= [hDtf , Hi := -ih—{hDt) . (5.64) 

a 

Remark 5.8. The commutator €2,3 gives the leading order term. This can be 
understood by a 'back-of-the-envelope' calculation replacing H by the leading terms 
from ijmil . I.e. H « Po + h^^^Pi + h^''^P2- We do not give this formal calculation 
here, since we do not justify this approximation. 

5.4.3. Control of ^^{Cixi o'^h\xi oVh) ■ The terms with derivatives in Ds are the 
most involved. 

Commutation vi^ith Tii. 

Since Xt commutes with W and Hi, we find 

[XtW, [xtW,Hi]] = XtlW, [W,Hi]] . (5.65) 

The inner commutator becomes 

[W,ni] Qi + Q2 , (5.66) 
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with 

Qi := {hD, - Ai)[W, a-^]ihD, - ii) , (5.67) 
Q2:=-{[W,Ai]a-\hD,-Ai) + {hDs-Ai)a-^[W,Ai]'j . (5.68) 

We calculate the double commutators separately 
[W,Qi] = (hDs ~ Ai)[W,[W,a'^]]ihDs - Ai) 

- { [W, ii] [W, ^hD, - ii) + {hDs - ii ) [W, a-2] [w, ii] } , 
[W,Q2] = -[W, [W,Ai]]a-^{hDs - ii) - [hD, - Ai)a'^[W, [W,Ai]] 

- [W,Ai][W,a-\hDs - ii) - {hDs - Ai)[W, a~^][W, Ai] 

+ 2[W,Ai]a-^[W,Ai] . (5.69) 

Remember that Ai{s,t) = — t(l — tK(s)/2), a{s,t) = 1 — tK{s). Therefore, 

[W, ii] = ^(^)]/2 = <2/ii/2-^o^ ^ [^^ [ly^ Ai]] fh'-^'02 , (5.70) 
[W, = t/ii/s-^Og ^ [w, [W, a-'^]] = t/i^^^e^^^ ^ (5 7^^ 

where, after a right multiplication by a cutoff function localizing in [0, to), the Oj's 
are bounded (pseudodifferential) operators commuting with the multiplication by 
functions tpit)- 

Using H5.70|) . H5.71|l . the Cauchy-Schwarz inequahty and that \t\ < 2/i^/^^* on 
suppxt, we find from H5.69|l . 

lixlo^h I [xtW,[xtW,ni]]xlo^h)\ 

<c[hi-'^-'\\{hD^~A,)xtxlo^kr 

+ h'-'^-''\\{hD^ - A,)xtxlo^h\\ WxtxioVkW + h'-''-*'\\xtxioVkr) 

< ch^/'-^^-''{{xm. I n{h)xtvh) + {h + h^/'-^')\\xm,r) ■ (5.72) 

Using the condition satisfied by S from (|5.7|l . and the localization estimate in s 
(|5.56|l . this implies (|5.58|) for the expectation value (xtfifh \ [xtW, [xtW,Tii]]x'ifi^h) ■ 
Commutation with H2- 

The commutation with Ti.2 is similar, but easier. 

[xtW, [xtW,n2]\ = -ihxl[W, [W, ^]{hD, - ii) - ^[l^,ii)]] 

= -ihxlW, ^]]{hDs - ii) + 2ihxl[W, ^][iy,ii)] 
+ ihxl^[W,[W,A^)]]. (5.73) 
Now, 9sa = — tK'(s), so the new terms to control are 

[ly, ^] = h'/2-H O5 , [ly, [ly, i?]] = h^-'H O, , (5.74) 
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with bounded^ O/s. Thus, for [xtW, [xtW,'H2]], we get 

KxloVh I [xtW,[xtW,n2]]xlo^h)\ 

< Ch^-^'i\\txtxiofh\\ \\{hD, - A,)xtxloV>h\\ + \\t\txloVh\\ \\tXtxloV>h\\) 

< C + ((;^,^^ I H{h)xtVh) + hWxtVhf) ■ (5.75) 

This implies for {xloVh \ [xtW, [xtW,'H2]]xlo^h) ■ 

Commutation with H3. 

Since W commutes with xt and Dt, we can calculate 

[XtW, [xtW,n3]] = h^W^ixu [xt,D^]] = -2h^W^\dtXt\^ 

= -2h'+^'w'\x[ijJ^)\'. (5.76) 

The expectation value of this term in the state oVh will be exponentially small, 
due to the normal Agmon estimates, H5.55|l . Explicitly, 

Kxlo^h I [xtW,[xtW,H,]]xloVh)\ 

= 2h'+''\{e"*/'^''\ioV, I w'e~^-'/''''\x[{j^^^re-'^'''\lo^,) 

<2/.i+2*||^^|p||e-2-/''^^^|x'i(^)nioo||e"*/''^'V/.f 

< Ch^+^^e-^"'''' = Ou„if(/i°°) . (5.77) 

In particular, (|5.58(l is satisfied for the term (xlfifh \ [xtW, [xtW,l-Lz\\X?s,a'-Ph) ■ 
Commutation with 7^4. 

When we calculate [xtW^, [Xi^^i ^4]]i we will use the discussion of the previous 
paragraph to conclude that if a derivative falls on xt , then the resulting expectation 
value becomes exponentially small. Thus, 



ixlo^h I [xtW,[xtW,ni]]xloVh) 

= ixloVh I - th[W, [W, ^]]xUhDt)xlo^H) + Ou„if(/i°") 



= -^Hxtxlo^h I - ih[W^ [W, J^]]{hDt)xtxlflVh) + Ounif(/i°°) . 

(5.78) 

From the formula — — , we see that all derivatives (in s) of — are uni- 

formly bounded on the support of xt- We therefore find 

[W, [W, ^]] = h'-^^O^ , (5.79) 
a 

where Oy is a bounded (pseudodifferential) operator. Thus 

I {xloVh I [xtW, [xtW, n^Wxlfl^h) I (5.80) 

< Ch^-^'{{xtvh I n(h)xm,) + hWxtVhW^) + Ou„if(/i°°) . 



''after multiplication by a cutoff function localizing in [0,to)i 
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This implies (15.58(1 for {xlfi^h \ [xtW,[xtW,TL4\]xlfl'Ph) , and therefore (|5.58|l is 
estabhshed for the expectation of Ci . 

5.4.4. Control of \{C2xlfiVh \ xloVh)- 

To estimate this term, we use similar calculations and the decomposition given in 
(|^:HT| and lEEU. 

Estimate of (C2Axl,ofh , xlfiVh)- 

For the first term VF[xt, [xt,H]]W^, we clearly get, using H5.55|l and as for l|5.77|l . 

{xlo^h\W[xu [XumWxlo^h) - -h'{xlo^h\W\dtXt?Wxlo^h) 

= Ounif . (5.81) 

Thus, H5.58|l holds for the expectation of €2,1- 
Estimate of {C2axlfi'Ph\xlflVh)- 

For this term, we notice that [VF, Tia] = 0, and calculate, using the Oj's introduced 
previously. 

C2,2 = Wxl[w,ni + n2 + 'Hi\ 

= Wxl({hDs - Ai)[W,a-\hDs - ii) - [W, Ai\a'^{hDs - ii) 

- {hDs - ii)a-2[M^, - ih[W, -^]{hD, - ii) 

+ ^h^[W,A,]-^h[W,^]{hD,)) 

= Wxl{{hDs - Ai)th^l'^-''0:i{hDs - ii) - t^h^'^"Oia-^{hDs - ii) 

- {hDs - Ai)a-H^h^/^-'Oi - ith^/^~''05{hDs - ii) 

_ ih^/^-H'^t^Oi ~ ih^/^-'OsihDt)) . (5.82) 
Here Os (which we used in the last line) and Og (that will be used below), are 

[W, = h^/^-'Os , W = th^^^-'Og . (5.83) 

a 

Therefore, we get (using the pseudo-differential calculus for showing that [Ds, Og] — 

\{xlo'Ph\C2axlo'Ph)\ 

< c{h^-^'h^-^'\\(hD, - A^)OgxU^h\\ WihD, - A,)xlo^h\\ 

+ h^-^^h^/^-''\\xl,VH\\ \\{hD, A,)xI,Vh\\ 

+ h'-'^h^/'-^'UhD, - A,)Ogxlo^h\\ llx^.o^'.ll 

+ h'-'^h'-''\\xloV,\\ \\{hD, - ii)4o^.|| + h'-''h'-''\\xlo^Hf 

+ h'-'^h'/'-'\\xlo^,\\\\hD,xlo^,\\) 

< Ch^'/^-^'diphlniph) + hWiphf) . (5.84) 
So, by (|5.7|l . (|5.58|) also holds for the expectation of €2,2- 
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Estimate of {C2,d.xlfifh\xlfifh)- 

For this last term the approach is equally direct. We calculate : 

^,n\ = [H?, {hD, - A^)]a-^{hD, - ii) + {hD, - A^)^,a-\hD, - ii) 
+ {hD,-Ai)ar^^,{hD,-A^)] 

- ^\{hDs - ii) - ih^[W, {hD, - ii)] 

-[*-,{hDtnW,K{s)] 

+ ih- [[W, Kis)], ^] (hDt) + ih^ [-, (hD,)] [W, k{s)] . (5.85) 

(X Ct Qi Qi 

The commutators in the above expression can be estimated ; 

[W,{hDs-Ai)]=-t[a-'[W.Kis)],hDs] +t[a-'[W,K{s)],Ai] 

= -iht-^ [W, k(s)] + hta-^[[W, k{s)], D^] 
-t^a-'[[W,K{s)],K{s)]/2 

= -ihf^[W,K{s)]+hta-^[W, [k{s),Ds]] 

-t^a-'[[W,K{s)],^{s)]/2 
= + th^/'^-^U(2) + t^h^-'^m^s) , (5.86) 

with bounded l^{j)'s when composed with t-cut-off functions. 
Thus, when we localize in {\t\ < Ch}''^-^}, we get 

[W^, {hD, - ii)] = h''-^-'Ui , (5.87) 
where Ui is uniformly bounded. Similarly, 

[W, a-'] = [W, = t\-' [[W, n{s)], ^^^'^ ~ ''^^'^ j 

= h}-^^h}-^'U2 , (5.88) 



(5.89) 



[-^,{hDt)]=ihd/-) . (5.90) 
[-, {hDtf] = ih{hDt)dt{-) + ihdt{-){hDt) = 2ih{hDt)dt{-) - h^d^i-) 

d Qi Cl CI Ct 

= {hDt)hUi + h^Us , (5.91) 



[[W,K{s)],^]=h'-^'Ue, (5.92) 
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Thus, we can estimate, after additional controls of commutators, 

\{xlflVh\C23xlfiVh)\ 

< cih'-'-^Wxlo^HW \\{hD, - A,)Wx^txloVh\\ 

+ h^-^'h^-^^\\{hDs - M)xl^^h\\ \\{hDs - A,)WxhloVh\\ 

+ e-'-^\\{hDs-A,)xloVh\\\\xlflVh\\ 

+ h^-^'h'-^^\\xloVH\\ \\{hD, - A,)WxlxlnVh\\ 

+ h^'-'-lxloVhf + h\{xlo^H I {hD,)U^[W, K{s)]xloVh)\ + h^'^-^WxloVuf 
+ h^/^-'h'-^^\\xl,M \\{hD,)WxWsflM+h'"^^\\xl,^uf) 

< C(h^'^''~^{^u I H^h) + h''^~'"^\\^uf 

+ h\{xlo^H\{hD,)Ui[W,n{s)]xl^^h)\) ■ (5.93) 

Remark 5.9. //, in (|5.93|l . we estimate [PV, k(s)] by the easy pseudodifferential 
result, i.e. use the bound 

\\[W,^{s)]\\^0{h^/^-^), 
we can continue H5.93|l to get 

\{xlo^h\C2^3xlo^h)\ < c{h'-H^,,\n^h) + h'-^yn\\') . 

This would only allow us to take e G (0,1/3) instead of e € (0,3/8) as claimed in 
Theorem \5.1\ In order to get the optimal range for e, we estimate the commutator 
[W, k{s)] .slightly better using the Agmon estimates' in s, i.e. H5.56|l . 

We write k{s) = fcmax — s^k(s), with k a smooth function with bounded deriva- 
tives of all orders. Now, 

[W, k{.s)] = [k{.s)s^,W] = k{s)[s, [s, W]] + 2k(s)[s, W]s + [k{s), W].s^ 

= Uih^'^' + U2h^/^-'s + Uah^^^-'s'' . (5.94) 

Therefore, using the estimate (I5.56|l . we get the inequality, 

\\[W,ii{s)]xlo^h\\ < Ch'^''''-^\\iph\\+Onniiih°°) , (5.95) 

and finally, H5.58|l follows (for the final term €2,3) by inserting H5.95|l in H5.93|l . 
This ends the proof of Proposition 15. 71 

6. Lower bounds in Grushin's approach 
In this section we finish the proof of Theorem ll.il 

Theorem 6.1. Let ^^'^\h) be the n-th eigenvalue ofH^h) and let zjxi\h) be the 
asymptotic sum given in (|3.4() . Then ^^"^^h) has z^\h) as asymptotic expansion. 

It is clear that Theorem 16.11 implies Theorem 11.11 From the estimates in the 
previous sections 0] and we know the low- lying eigenfunctions to be localized in 
phase space, around the following set (where we denote by ^ the variable in phase 
space dual to s) : 

• s e (-2/zi/8-Ei ^ 2/ii/8-^i) ; 
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. e e (-2/1^-^-1/2 +^q/,-1/2 ^ 2/l^^'"l/2 +^g/j-l/2) 

We will choose ci = £2 ?? , £3 = 3/8 — 77 , where 77 < 1/8 . Thus, in the remainder 
of the section, 77 will be a fixed (small) constant satisfying, 

77 6(0,1/8). (6.1) 

Remark 6.2. 

We only need a space localization in t. There is no localization in the correspond- 
ing frequency. This is pleasant, since it avoids possible complications due to the 
boundary condition. 

More precisely, the localization of is analyzed in the following lemma. 
Lemma 6.3. 

Let M > 0, Hq > and let x be a standard cut-off function: 

X&C^m, x{x)^l, onanbd. of[-^,+^], suppx C (-2, 2) . (6.2) 

Then for all K > there exists Bk > such that if {uh)heio.ha) is a family of 
normalized eigenf unctions of TL with eigenvalue fi{h) satisfying 

fi{h) < Ooh - Cifcmax/i^/^ + , (6.3) 

then, with rj from (|6.1|l . 

Proof of Lemma \6.Sl 

Define X2 = 1 ^ X- From Theorems 14 . 1 1 and 14 . 91 we know that 
So it suffices to prove that 

= Ounif(/i°°) . (6.5) 

By writing 

= v(^)x(^)|x( "^"' ) + V2( "^°' )|X( — ).; 
and appealing to Theorem lS.ll we get (|6.5I) and thereby Lemma [6.31 □ 



2 



Proof of Theorem \6.1\ 

For a definite choice of x (fixed once and for all) as in Lemma 16.31 and Uf^ an 
eigenfunction of Ti, satisfying H6.3|) . define 

^'^ = xij^)xij^)x C ,3/8-. 
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Calculations will from now on be carried out in the variables (cr, r) . All functions 
considered will be localized on a scale of order h^^ in the (cr, T)-variables. This 
implies (in particular) that they are localized to a tubular neighborhood of size 
/ji/2-'j near the boundary in the original coordinates x £ fl. The natural measure 
in (cr, r)- variables, inherited from L^{rt,dx) by implementing unitarily the change 
of coordinates, is {1 — h^^'^TK{h~^^^a))dadT. However, due to the localization of our 
functions (and the boundedness of k) we can replace this measure by dadr (since 
h^^'^Tfi{h^-^/^a) — 0{h^^^^''^) on |r| < Ch"^^), without changing our estimates up 
to multiplicative /i-independent constants. Therefore we can (and will !) do all our 
estimates by choosing the norms in i^(R x M+, dadr) or in £(L^(R x R+, dadr)). 
Thus all i^-norms below refer to L^(MxM_|_, dadr) and similarly for operator norms. 

With these conventions we have, using Theorems 14. II and I^TTl for all eigen- 
functions Uh corresponding to eigenvalues ^i{h) satisfying (|6.3I) . that the correspond- 
ing if^h (given in IjG.GIl ) satisfies (with error terms Ounif (^°°) uniform for eigenfunc- 
tions Uh as long as H6.3|l is satisfied), 

llV^ftlU^ = l + anif(/i°°) , and L4,h=^h+0^^,i{h^) , (6.7) 

for all 

L = xij^Mj^mh^D^) , (6.8) 

with X satisfying (|6.2I) and with rj from H6.1|) . Let us fix an Lq as in (|6.8|) in the 
rest of this Section. 

Let i/harm be the harmonic oscillator on L^(M) defined by 

iJharm 3Ci Ve~oDl + Ci ^ , (6.9) 
(compare with Section |3J|. Clearly, -ffharm has eigenvalues 

e, :=Ciey'yif (2^+1) , (6.10) 

with £ £ N. Let V{ be the corresponding (unique up to scalar multiple) normalized 
eigenfunction. For N £ N, the value CiB^'^VGfci^A^ is right in the middle be- 
tween two eigenvalues {cn-i and cat). We define the vector space Vn C L^(M) as 
the space spanned by eigenfunctions of H\is,rn\ corresponding to eigenvalues below 
ClQy^^/W2N, i.e. 

Vn ■■= Ranljg^^^g,i/4^/g^^](-ffharm) = Spanjwo, • ■ • ,-^^-1} • (6.11) 

Clearly, dim Vn = N. 

Similarly, we define UN{h) C L'^{il) as the spectral subspace attached to the 
interval /Ar(/i), with 

iNih) = ( - C30, Go/l - fcmaxCi/l3/2 + CiOl^^ ^Gk^N . (6.12) 

Let '■ i^(R) Vn and : L'^{Q) UnQt) be the orthogonal projections. 
We define a linear map M[^\h) from Vat to UN{h) by 

(where was defined in (|8.2(l ) and extended by linearity. The number Mq is 
chosen fixed, but suffficiently large — the choice Mq = 10 would suffice. 
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Furthermore, we define a linear map Mf\h) from UN{h) to Vn by 

Mf\K)uh = Uvj^RoLoijjh , 
wfiere iph is defined from Uh by H6.6|l . We will prove the following lemma. 
Lemma 6.4. Let N eN \ {0}. Then there exists Hq > 0, and rj > (as in (|6.1|) 
and used in (|6.t)|) ) such that for all h < Hq, A4\^\h) and M2^\h) are bijective. 

We will prove Lemma 16.41 below. First we apply it to finish the proof of Theo- 
rem 

Lemma |6.4I implies that, for sufficiently small h, 

dimllNih) = dim Fat N . 
But Corollary 1^31 describes N distinct points of Spec(7i) below 

This finishes the proof of Theorem 16. II □ 
Proof of Lemma |6'.^[ 

We only need to prove that M[^\h) and M2^\h) are both injective. Injectivity 
of A4[^\h) is clear from Sectional so we only consider injectivity of ^A2^\h). 
The key to the proof of injectivity of A42^\h) is the following lemma. 

Lemma 6.5. There exists 770 > such that if rj < rjQ in (|6.1|l . then there exists 
C > such that for all normalized eigenf unctions Ufi G Uj^{K) with corresponding 
eigenvalue ^{h), we have 

-i7harm)i?^i0^/.|| < Ch^'^'' , (6.13) 

where iph is related to by (|6.6(l and ^{h) is defined by 

y[h) := h-"\^J,{h) - {Boh - k^^^Cih^/^)} . (6.14) 
Proof of Lemma \E~^ With P from H3.1()(l and 

A(/i) = /i-i(/^(/i)-eo/i) , (6-15) 
we have (using Theorems 14 . 1 1 IT!^ and l5.1|l . uniformly for normalized eigenfunctions 
Uh € UN{h), 

(P ~ Xih))^h = Onniiih^) , and (P - A(/i))L^ft = Ou„if (/i°") . (6.16) 

In the rest of the proof of Lemma 16.51 we will often have estimates like H6.16|l . 
We will generally not repeat the phrase 'uniformly for normalized eigenfunctions 
Uh € Uj\j{hy, but the estimates are meant to have such uniformity. 

Using (|6.16|) and the notation from (|3.24|) . we get: 

'P ~ \{h) (L^^h\ 

Bo J\ ^ J yRoLoi'hy 



+ Ounif(/l°°) . (6.17) 



Furthermore, with £q from (|3.24|) . 



(6.18) 



Here we have introduced 

dP={P- Po) - X{h) . (6.19) 
In order to proceed we need a bound on the matrix in {•} in (|6.18|l . 
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Lemma 6.6. 

There exists a constant C > such that ( with rj from H6.1(l ) 

\\E^{P - Po)io|| + ||i?o - PM\ < C/i^/s-" . (6.20) 

More precisely, with Pi, P2 and P3 from (|3.16|l - ()3.18|l . there exists Nq E N such 
that, for rj satisfying (|6.1|) . 

\\Eo{P~Po)Lo - Eoih^^^'Pi + h^'^P^ + h^'^P,)L,\\ < Chy^-^«^ , 

\\Ro (P - Po)Lo - Ro {h^^'Pi + h'/^P2 + h^^*P3)Lo\\ < Ch'/'-''°'^ . (6.21) 

Furthermore, 

h'/'\\EoPiLo\\+h'/'\\RoPiLo\\ < Chi-^^ , 

h'/^EoP2Lo\\+h'^^R-P2Lo\\ < Ch'/'-^^ , 

h^'/^WEoPsLaW + /i^/*||i?o"P3io|| < C;i3/4-5, _ 22) 



Proof of Lemma \6.fA 

With a and 0,2 from H3.9fl and omitting the tilda's on the a's, we have by definition 

P = a-' [{t + Co) + h^^^D, - r(l - 02)] a"' [(r + Co) + - r(l - 03)] 

+ a-^DroDr - 60 , (6.23) 
and therefore 

P-Po = [a-^ [{t + Co) + h^'^D„ - r(l - 02)] a'' [(r + ^o) + h^^^D, 

- r(l - 02)] -{r + ^0)'} - ^Dr^ + dr{^) . (6.24) 

We win use the property that L locahzes to {r < 2h^^, \a\ < 2h^^} and that 
Eg , EqDt , R^ , and R^ arc bounded. We introduce f '■= {t + ^q) — t{1 — 02) 
and calculate 

a-\f + h^/^D,)a-\f + h^'^D,) - (r + 

= { f/a^ - (r + Co)') + ih^/'^f + 2a- V/^'/'^. + ^^h^"D, + a-2/^3/4i?^ ^ 
Thus, 



P - = Qi + Q2 + Q3 , 



where 



Qi - /V«^ - (r + Co)^ + ^h'/'^f + dr (^) , 

a"* a 

Q2 = -iDr , 

a 

Now, on {r < 2/1"'', |cr| < 2ft.-''}, we have 

a= l + 0(fti/2-';) ^ 1 - 02 = ©(ft'/'-") . (6.25) 



MAGNETIC BOTTLES 



43 



Therefore, 

dra = 0{h^'^), dr(—)=0{h). 

a 

Thus, Qi = ©(/ji/^-^'f), so 

\\EoQiL4 + \\R^QM=0{h^-^^) . 
Furthermore, for all j > 0, 

Dix{h^a)x{h^D,) = x{2h^<j)Dix{h^a)x{h^D,) , (6.26) 

and 

\\Dixih'^a)x{h^D,)\\ < c,ih~y . (6.27) 
Using and if^TTjl . we get 

||i^oQ2io|| + ||i?o"Q2io|| =0(/l'/'), 



This finishes the proof of 

The more precise estimates, (|6.21() and (|6.22|) . follow in the same manner, using 
that on {t < 2h-'^, \a\ < 2/1^''}, 

a = 1 - h^/^r{K{0) - ^h^/^a^n"{0) + ©(/i^d/s--)))) ^ 

a2 = 1 - /i^/Vi (k(0) - ^h'/*a^K"{0) + ©(/i^d/s-"))) , 

instead of H6.25|l . We omit the details. 

This finishes the proof of Lemma 1^?^ □ 

Combining (|6.17(l and (|6.18() . we get : 

Lo^h = RtRoLo^h - Eo{dP)Loi;h + Ou„if(/i°°) , (6.28) 

and 

R^{dP)Loi,h = Ou„if(/j°°) . (6.29) 

We now introduce an additional localization through an operator L as in H6.8|l . 
which is chosen (slightly 'bigger' than L^, i.e.) such that 

LLo^ Lo + Oih°°). (6.30) 

We observe that H6.29|l is also valid with Lq replaced by LLq and, applying the 
(uniformly) bounded operator R^{dP)L to (|6.28|) . we obtain 

R^{dP)LR+R^Loil>h - RoidP)LEo{dP)Loi'h = Ou„if(/i°") • (6.31) 
We again apply Lemma 16.61 (|6.31|) and the comparison estimates (|6.34f) and 



(to be proved below) and obtain that, for all 5 > 0, there exists 770 £ (0, i), 
such that if 77 < 770 , then 

- R^{h^'^Pi + h^'^P2 + /i^/^Pa - X{h))LR+R^Loi;h 

+ Ro{h^/^Pi)LEo{h^/^Pi)R^RoLoi^h 
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Using H6.29|l . the rapid decay of the function uq, the support properties (in r) of 
L and the pseudo-differential calculus in the a variable for controlling commutators 
(in order to push L to the right), we finally get : 

- i?(7(/i'/*'Pi + h^^^P2 + h^/^P3 - X{h))R+RoLoiJh 

+ Ro{h''/''Pi)Eo{h''/^Pi)R^RoLo^h 

= Oun-Ahn + O^n-Ah'^'-'m^R^LoM ■ (6.32) 

We get H6.13|l from (|fi.32|l by calculations similar to those leading to the expressions 
for El, E2, E3 in Subsection 13.41 We just recall that 

RqPiR^ = , -RqP^R^ - k(0)Ci , 

and 

— Rq P3R0 + Ro Pi^oPiRq — —Hha-cm ■ 

This finishes the proof of Lemma 1^31 □ 

We now compare (as already used above) various norms and observe : 
Lemma 6.7. 

Let G N. There exists c > and ho > such that ifiph is associated (as in (|6.6(l ) 
to a normalized eigenfunction Uh of TL with Uh € UM{h), then for all h £ (0, /iq], 

||i?^Lo^,,|| - < ll^^ll < WRoLoMl + ch'/^ . (6.33) 

Proof. 

Since clearly \\Rq\\ = 1, we get from (|6.7|l . 

This implies the first inequality in Lemma 16.71 To get the second inequality, we 
apply H6.28|l . Lemma RTHl and (|6.7(l . and get 

= R^RoLo^h + Ounif (/i'/'^") . (6.34) 
Since \\R^\\ — 1 and tph satisfies (|6.7f) . this implies the lemma. □ 

Using the self-adjointness of the harmonic oscillators, we get the following propo- 
sition. 

Proposition 6.8. Let N E N. There exist ho > and C > such that if 
{ti'{h))he{o,ho] eigenvalue of Ti satisfying with /i(ft.) G lN{h) (see (jfi.f 2|l ). 

then v{h) (defined by (16.14(1 ) satisfies 

y{h) e U.^-^ea + [-C/i^/ie^ +Ch^/^''] . (6.35) 

Proof. Using Lemma l6 . 71 above. Lemma 16.51 implies that 

dist(z.(/i),Spec{i/ha„„}) = Ou„if(/i'/'') . 

□ 

Lemma 6.9. Let iV e N. There exists ho > such that if h £ {0,ho], then 
dimUNih) = N. 
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Proof. Wc know from Section 13 that dimUN{h) > N. In order to prove Lemma 
we only have to prove that the eigenspace attached to some interval I'ih) G [ee — 
Ch^^^^,ei + C/i^/^^], with from (|6.10() . and I < N , is necessarily of dimension 
< 1. If it was not the case, let Mi.^, U2.h be normalized orthogonal eigenfunctions 
corresponding to eigenvalues fJ.i{h) and H2{h) in the interval 

Qoh - Cik^^^h"/^ + h"^[et - e, + C/i^/i^] , 

for some £ — l{h). Let V'l./u "^^.h be defined as in (jfi.fifl and let i^i, V2 be as in H().14|l . 
Let ei and vi be as in (|6.1U|) and below. 

For a,b G C {a,b will depend on h) with |ap + |5p — 1. We have, using the 
almost orthonormality of t/^i.h, 4'2.h and (|6.7|l . 

1 + Ouniiihn = \\atPi,h + b^2.h\\^ < \\RoLo{atlJi,h + b^2,hW + ©unif (/i'/'"") • 

(6.36) 

With £ = £{h) as above, we may choose a, h such that 



u^((T)i?o Loiaipi^h + bip2,h) da = . (6.37) 

X) 

Lemma 16.51 implies that 

{a - Hharm)Ro Loiaipl^h + b'lp2,h) 
= aifiih) - ffharm)i?o"ioV'l./i + b{iy2{h) - i?harm)i?(7ioV'2,h + ©unif (^l^^^^ ) 

= Ou„if(/i'/'') . (6.38) 
Using (|07|) . imphes that 

||i?^Lo(aV-'ij, + 6V2,/«)|| = Ounif(/^'/''^) , (6.39) 

which is in contradiction to (|6.36|) . This finishes the proof of Lemma [6.91 □ 

Thus, for sufficiently small h, UN{h) ~ Spa,ii{u\p}^^-^^. 
Lemma 6.10. Let N £ J\f. There exists ho > such that 

for all h < ho and all j £ {1, ... , N}. 

Proof. By induction it suffices to prove the lemma for j = N. By Lemma f6. 51 and 
the spectral theorem there exists £{h) e {0, ... — 1} such that (with ip^ being 
associated to ulj^' as in H6.6|) 'l 

Suppose £{hn) < iV — 1 for a sequence {hn} with hn \ 0. Then u\l^'' e UN-i{hn) 
and therefore dim ?7jv-i(/in) > N, in contradiction to Lemma 16.91 Thus, 

RqLoiJjn - {VN-I, Rq Lq'4jn)vn~i = Ounil{h^^^^) , 

and therefore 

X^uW = {vN^l,RoLoMvN-l+0^nit{h^^^'') ■ 

Lemma l6 . 1 01 now follows from Lemma 16.71 □ 

The injectivity of M.2^'' clearly follows from Lemma [6. 101 This finishes the proof 
of Lemma 16.41 □ 
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Appendix A. On an important family of ordinary differential 

EQUATIONS 

Let us recall for the comfort of the reader the main properties (mainly due 
to |DaHe| and PeHt] ). concerning the Neumann realization of H^-^ in L^(]R+) 
associated to + + We denote by //^^•'(C) the lowest eigenvalue of H^'^ and 
by (p(^ the corresponding strictly positive normalized eigenfunction. More simply 
we will write /i(^) instead of A^^HO this appendix. It has been proved that 
the infimum inf^gR inf Spec(iJ^'^) is actually a minimum. Then one can show that 
there exists < such that /i(^) decays monotonically to a minimum value Oq < 1 
and then increases monotonically again. So it can be proved that : 

Go = inf (inf Spec(i7^'«)) = inf Spec{H^'^°) , (A.l) 
and moreover that: 

Go = Co ■ (A.2) 

It is indeed proved in [ DaHej that 

M'(O = [M^)-e>«(0f . (A.3) 

From (|A.3p . we get that 

m"(6) = -26¥'L(0)>0. (A.4) 

We will write uq instead of , and define the constant Ci by . 

Let us now recall some formulas appearing in |BeStj . Define Mk to be the fc'th 
moment, centered at —£,o, of the measure Uq{x) dx : 

Mk = ix + ^o)''ul{x) dx . (A.5) 



These moments were calculated in jBeSt| . 
Lemma A.l. 

The first moments can be expressed by the following formulas : 

Mo = l, Mi=0, ^^2-^, M3 = ^>0. (A.6) 

We will also need a few other results on the model operator. 

Proposition A.2. 

We have the following identities 

j^[2r{r + Co)^ - r\r + ^,)]ul{r) dr = ^ ^ ^ > , 
■ r f \n ( \A "o(0) 3 

I I Uo[T)DrUo[T)dT = = --^Cl . 

Proof. 

The first identity clearly follows from the known moments of Uq and The 
second identity follows from partial integration. □ 

Proposition A.3. 

For z e M, let E{z) be defined as the ground state energy of the Neumann realization 
of 

H{z) = -^ + {T + i, + z)\ 
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on L^(M+). Then E[z) is a smooth function and satisfies 

E"{0) = 2(1 - 4/2) , (A.7) 

with I2 from (|3.37|l . 
Furthermore, 

E"{Q) = 6Ci^/e^ . (A.8) 

Proof. 

By analytic perturbation theory, E{z) is analytic and there exists an analytic func- 
tion R 9 z (j){z) e such that 

U{z)\\ = 1 , H{z)cj){z) = E{z)cj){z) , </.(0) = . (A.9) 

By differentiating the identity = 1 twice with respect to z, we find 

25ft(0'(O) I lio) - , -||0'(O)||2 = 5R(0"(O) l^/o) . (A.IO) 

From the equation H{z)(j){z) ~ E{z)(f){z), and the fact that E{z) is minimal at 
z = 0, we get, with Pq = H{0) - Gq , 

Po(/)'(0) = -2(t + $oH, 
which imphcs with P^^ from (|A.13|I and (|A.14p . since uq -L (t + Co)wo (by HAjp ). 

0'(O) - -2Po-i((T + eo)"o) , (A.ll) 

for some c e iM. Finally, differentiating the relation E{z) — {(j){z) \ H{z)4){z)) twice 
gives us the formula : 

^"(0) = 2eo5R(0"(O) I uo) + 83?(</)'(0) I (r + Co)"o) 

+ 2((/.'(0) |ff(O)0'(O))+2. (A.12) 

Upon inserting HA.10|I and (|A.11|I in (|A.12|I . we get (|A.7p . The final identity, HA.8|) 
is a rephrasing of (|A.4|I . □ 

We also have the following easy observation : 

Lemma A. 4. 

Let he the operator from Suppose 4> £ 5(R), then R^4> G 5(M x M^) • 

Proo/. 

This is an easy consequence of the regularity and decay of mq . □ 

Finally, we will need the following mapping properties of the regularized resol- 
vent. 

Lemma A. 5. 

Let Pq be the Neumann realization of 

~^ + ir + Cof-eo, 

on L'^(R^). For ^ _L uq we can define Pq^4> as the unique solution f to 

Pof^<f>, f±uo. (A.13) 

Let Pq'^ £ C{L'^{R+)) be the regularized resolvent : 

""^'^A'p-^s 1t"° (^-^^^ 
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(and extended by linearity). Then Pq^ is continuous from iS(M_|_) into »S(M+) . 

Moreover, for any a > 0, Pq^ is continuous in Z/^(M+ ; exp— ar) . 

Proof. 

Using the local regularity up to the boundary of Pq, one first gets that Pq^ sends 
<S(R+) into C°°(M+). For the control at oo, one then observes, after cutting away 
from 0, that the problem is reduced to the analysis of inverting the harmonic 
oscillator — ^ + {t}'^ — 6o on iS(Ili), which is a standard result. 
For the last statement, we can also observe that, for any real a, the operator 

exp-aVl + t2 • ( — f- + (r + ^of - 6o)"^ -expavT^Tr^ 

extends continuously on L'^{R) and S{R). 
One can also show by the same technique that 

t^u''o\t) e L'^(R+ ; exp -ar) , for all a > , and for all j, k . (A.15) 
With additional work, one could actually get a better decay. □ 

Appendix B. Coordinates near the boundary. 

It is convenient in most calculations to straighten out the boundary by a co- 
ordinate transformation. This, quite standard, procedure, will be defined below. 
Let zq G on and let i be the length of the boundary dVl and / — |, |] . Let 
M G C°°{I;dn) be a parametrization of such that M(0) = zo and s is the 
distance inside dCl between M(s) and zq. We denote by 

T{s) := M'{s) , 

the unit tangent vector of dSl at M{s) and the scalar curvature by k{s), which can 
be defined by 

T'{s)=k{s)u{s) , 

where f(s) is the interior normal unit vector of 90 at M{s). 
Moreover the parametrization is chosen positive : 

det (T(s),i/(s)) = 1, Vs e / . 

For any z G O , we denote by t{z) the standard distance of z to dfl : 

t{z) = inf 1^; — w| . 

So, there exists to > and a diffeomorphism of class C°° : 

^ 51/(2,) x(0,io), (B.l) 

such that iIj{z) = w = {s{z),t{z)) and \z — M{s{z)) \ — t{z) . 

We have denoted, for small enough e, by fig the tubular neighborhood of dO. : 

9., := {z e n; dist(z,aO) < e} 

and is the circle of radius r is identified with [— 7rr, 7rr[ . 
So we have the identity 

Z = M{s{z))+t{z)u{s{z)), e Oeo • (B.2) 

From this equality, it is easy to check that 

T{s{z)) = [1 - t{z)K{s{z))]Vs{z) and iy{s{z)) = m{z) . (B.3) 
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So for all u G H^{Q) such that supp(u) C il-sg , 
f \{hD^ ^ A)u\^ dz = 

J UJ 

\{hDt - A2)v\'^ + (1 - tn{s))-^\{hDs - Ai)v\'^] (1 - tK{s)) dw (B.4) 



K 



and 

' \u\'^dz^ [ , (B.5) 



J K 

with = u(ip~^{w)) , K = I x [O^to) , w — (s, <) and dw — dsdt . 

The magnetic potential A satisfies 

Ai ds + A2dt^ Ax dx + A2dy . 

So 

dA OA 

[-^(w) - -^H] ds/\dt = B{z) dx /\dy = B{w)[l - tK{s)] dshdt , (B.6) 
with ^j{z) = w and B defined as : 

B{w) = B{z) . (B.7) 

This gives : 

- = Bi^^'i^W - Hs)] = Bit, - Hs)) ■ (B.8) 

Then we get the identity between differential operators 

(hD, - A)^ = a-%hDs - M)a-\hD, - Ii) + {hDt ~ A2)a{hDt ~ I2)] , (B.9) 
where a{w) = 1 — tK{s). 

The usual Hilbert space L'^ifltg) is transformed to L'^{K;adw) . 

In the new coordinates and using a gauge transform, we can always assume that 

the magnetic potential has no normal component in a neighborhood of dfl : 

I2 = . (B.IO) 

In this case, we have : 

dtAi^~B{t,s)il-t,i{s)) , (B.ll) 
where B was introduced in (IB. 711. 
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